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[i]  The  Cryogenic  Infrared  Spectrometers  and  Telescopes  for  the  Atmosphere  (CRISTA) 
instrument  measured  stratospheric  temperatures  and  trace  species  concentrations  with 
high  precision  and  spatial  resolution  during  two  missions.  The  measuring  technique  is 
infrared  limb-sounding  of  optically  thin  emissions.  In  a  general  approach,  we  investigate  the 
applicability  of  the  technique  to  measure  gravity  waves  (GWs)  in  the  retrieved  temperature 
data.  It  is  shown  that  GWs  with  wavelengths  of  the  order  of  100-200  km  horizontally  can 
be  detected.  The  results  are  applicable  to  any  instrument  using  the  same  technique.  We 
discuss  additional  constraints  inherent  to  the  CRISTA  instrument.  The  vertical  field  of  view 
and  the  influence  of  the  sampling  and  retrieval  imply  that  waves  with  vertical  wavelengths 
~3-5  km  or  larger  can  be  retrieved.  Global  distributions  of  GW  fluctuations  were 
extracted  from  temperature  data  measured  by  CRISTA  using  Maximum  Entropy  Method 
(MEM)  and  Harmonic  Analysis  (HA),  yielding  height  profiles  of  vertical  wavelength  and 
peak  amplitude  for  fluctuations  in  each  scanned  profile.  The  method  is  discussed  and 
compared  to  Fourier  transform  analyses  and  standard  deviations.  Analysis  of  data  from  the 
first  mission  reveals  large  GW  amplitudes  in  the  stratosphere  over  southernmost  South 
America.  These  waves  obey  the  dispersion  relation  for  linear  two-dimensional  mountain 
waves  (MWs).  The  horizontal  structure  on  6  November  1994  is  compared  to  temperature 
fields  calculated  by  the  Pennsylvania  State  University  (PSU)/National  Center  for 
Atmospheric  Research  (NCAR)  mesoscale  model  (MM5).  It  is  demonstrated  that  precise 
knowledge  of  the  instrument’s  sensitivity  is  essential.  Particularly  good  agreement  is 
found  at  the  southern  tip  of  South  America  where  the  MM5  accurately  reproduces  the 
amplitudes  and  phases  of  a  large-scale  wave  with  400  km  horizontal  wavelength.  Targeted 
ray-tracing  simulations  allow  us  to  interpret  some  of  the  observed  wave  features.  A 
companion  paper  will  discuss  MWs  on  a  global  scale  and  estimates  the  fraction  that  MWs 
contribute  to  the  total  GW  energy  (Preusse  et  ah,  in  preparation,  2002).  IndexTerms:  3334 
Meteorology  and  Atmospheric  Dynamics:  Middle  atmosphere  dynamics  (0341,  0342);  3360  Meteorology 
and  Atmospheric  Dynamics:  Remote  sensing;  3384  Meteorology  and  Atmospheric  Dynamics:  Waves  and  tides 
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1.  Introduction 

[2]  Gravity  waves  (GWs)  are  important  for  the  middle 
atmosphere  dynamics  and  composition.  Originating  from 
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lower  altitudes,  GWs  propagate  upwards  into  the  strato¬ 
sphere  and  mesosphere  where  they  deposit  momentum  and 
energy  when  breaking.  Therefore  GWs  strongly  influence 
the  global  circulation  and  are  thought  to  be  a  major  cause  of 
global-scale  oscillations  such  as  the  Quasi  Biennial  Oscil¬ 
lation  (QBO)  [Dunkerton,  1997,  and  references  therein]  and 
the  Semiannual  Oscillation  (SAO)  [e.g.,  Mayr  et  ai,  1998]. 

[3]  However,  most  current  global  circulation  models  are 
not  able  to  explicitly  resolve  GWs  because  of  the  coarse 
spatial  resolution  of  the  models.  Even  models  which  can  be 
run  for  case  studies  with  sufficient  resolution  to  resolve 
GWs  have  to  parameterize  GWs  for  long-term  runs  due  to 
computing  constraints.  The  importance  of  a  good  represen¬ 
tation  of  GWs  to  obtain  realistic  global  circulation,  mean 
wind  and  temperature  structures  as  well  as  QBO  and  SAO 
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patterns  has  been  highlighted  by  many  authors  [e.g.,  Boville, 
1995;  Manzini  et  al,  1997;  Manzini  and  McFarlane,  1998; 
McLandress,  1998]. 

[4]  A  fundamental  question  in  this  context  is  what  and 
where  the  major  wave  sources  are.  Waves  excited  by  flow 
over  mountains  (mountain  waves,  MWs),  have  zero  phase 
speeds  relative  to  the  ground  (at  least  as  long  as  the  wind 
speed  does  not  change  strongly  with  time)  in  contrast  to  GWs 
from  other  sources  such  as  convection.  Because  nonoro- 
graphic  GWs  have  different  impacts  on  the  synoptic-scale 
dynamics,  they  are  commonly  treated  by  different  parameter¬ 
ization  schemes  [Boville,  1995;  Manzini  and  McFarlane, 
1998].  Thus  there  is  interest  in  determining  the  sources  of 
observed  GWs  and  in  estimating  the  fraction  that  orographic 
GWs  contribute  to  the  total  GW  flux  entering  the  middle 
atmosphere  on  global  scale  [Nastrom  and  Fritts,  1992]. 

[5]  Breaking  GWs  not  only  deposit  momentum  but  also 
induce  turbulence  and  thus  enhance  turbulent  mixing  in  the 
atmosphere.  This  influences  the  chemical  composition  espe¬ 
cially  in  the  upper  mesosphere  where  wave  breaking  is 
intense  [Garcia  and  Solomon,  1985;  Summers  et  al.,  1997]. 

[6]  Furthermore  GWs  might  play  an  important  role  in  the 
ozone  depletion  in  the  Arctic  winter  and  spring.  While 
temperatures  in  the  Arctic  on  global  scales  normally  remain 
too  warm  to  form  polar  stratospheric  clouds  (PSC)  and 
thereby  activate  chlorine  radicals,  temperatures  may  drop 
locally  below  the  frost  point  in  the  presence  of  strong  MWs 
as  observed  in  the  lee  of  the  Scandinavian  mountains 
[Dornbrack  et  al.,  1999].  Given  the  high  number  of 
mountain  ridges  at  high  latitudes  in  the  northern  hemi¬ 
sphere,  such  events  should  occur  rather  frequently  and 
therefore  be  a  significant  source  of  PSCs  on  synoptic  scales 
[Carslaw  et  al.,  1999a,  1999b;  Dornbrack  et  al.,  2001; 
Dornbrack  and  Leutbecher,  2001]. 

[7]  In  recent  decades  most  observational  studies  of  GWs 
have  used  ground-based  data  from  a  limited  number  of 
locations  [e.g.,  Fritts,  1984;  Tsuda  et  al.,  1991;  Eckermann 
etal,  1995',  Allen  and  Vincent,  1995].  Even  large  coordinated 
networks  of  ground-based  and  in  situ  measurements,  such  as 
the  Planetary  Scale  Mesopause  Observing  System  (PSMOS) 
and  the  Stratospheric  Processes  And  their  Role  in  Climate 
(SPARC)  programs,  still  cannot  provide  true  global  coverage. 
Global  coverage  is  now  provided  by  a  rapidly  increasing 
number  of  GW  studies  from  satellites.  To  date,  GW  signals 
have  been  analyzed  on  global  scales  using  temperature  data 
acquired  by  the  Limb  Infrared  Monitor  of  the  Stratosphere 
(LIMS)  [Fetzer,  1990;  Fetzer  and  Gille,  1994, 1996],  63  GHz 
limb  radiances  measured  by  the  Microwave  Limb  Sounder 
(MLS)  [Wu  and  Waters,  1996a,  1996b,  1997;  McLandress 
et  al.,  2000;  Jiang  and  Wu,  2001],  temperature  data  from  the 
Cryogenic  Infrared  Spectrometers  and  Telescopes  for  the 
Atmosphere  (CRISTA)  [Eckermann  and  Preusse,  1999; 
Preusse  et  al.,  1999,  2000],  and  temperature  retrievals  from 
the  meteorological  program  for  the  Global  Positioning  Sys¬ 
tem  (GPS/MET)  [Tsuda  et  al,  2000;  Nastrom  et  al,  2000]. 
In  addition,  in  the  upper  mesosphere/lower  thermosphere 
data  from  the  Wind  Imaging  Interferometer  (WINDII)  were 
investigated  for  GWs  [Wang  et  al.,  2000a,  2000b,  2001a, 
2001b].  GWs  have  also  been  investigated  in  saturated  radi¬ 
ances  taken  by  the  SPIRIT-3  radiometer  on  board  the  Mid¬ 
course  Space  Experiment  (MSX)  satellite  for  two  case  studies 
[Dewan  et  al.,  1998;  Picard  et  al,  1998]. 


[s]  In  case  of  the  LIMS  and  the  CRISTA  instrament  GW 
data  are  inferred  from  optically  thin  emission  limb  measure¬ 
ments.  Other  instruments  using  the  same  technique  are  the 
Cryogenic  Limb  Array  Etalon  Spectrometer  (CLAES) 
[Roche  et  al,  1993]  and  the  SABER  (Sounding  of  the 
Atmosphere  using  Broadband  Emission  Radiometry)  instru¬ 
ment  [Mertens  et  al,  2001].  In  addition,  the  upcoming 
mission  of  the  High-Resolution  Dynamics  Limb  Sounder 
(HIRDLS)  on  board  the  Earth  Observing  System  Chem 
(EOS-Chem)  satellite  will  provide  similar  global  measure¬ 
ments  at  a  rate  of  ^8000  profiles  a  day  with  global  coverage 
[Barnett  et  al,  1998].  The  high  spatial  resolution  of  these 
data  ought  to  make  this  instrament  very  well  suited  for  GW 
studies. 

[9]  This  paper  provides  a  detailed  assessment  of  GW 
fluctuations  (and  MWs  in  particular)  observed  in  temper¬ 
ature  data  acquired  by  CRISTA  infrared  limb  scan  obser¬ 
vations.  Many  of  the  results,  however,  are  general  and 
should  apply  to  other  instruments,  such  as  HIRDLS.  To 
interpret  GW  patterns  in  data  from  such  instruments,  a 
detailed  knowledge  of  the  sensitivity  to  different  horizontal 
and  vertical  wavelengths  is  essential.  In  section  2  we  derive 
an  analytical  approximation  for  the  sensitivity  of  CRISTA  to 
GWs  and  compare  this  to  numerical  results  derived  from 
high  spatial  resolution  radiative  transfer  calculations.  We 
demonstrate  that  a  two-dimensional  consideration  of  ray 
paths  through  the  GW  stractures  is  imperative.  The  results 
prove  that  CRISTA  can  detect  waves  with  horizontal  wave¬ 
length  larger  than  100-200  km. 

[10]  In  section  3  we  describe  the  CRISTA  experiment,  data 
acquisition  and  analysis.  The  analysis  of  GWs  in  acquired 
temperatures  employs  a  combination  of  Maximum  Entropy 
Method  (MEM)  and  Harmonic  Analysis  (HA),  which  is  well 
suited  for  the  analysis  of  monochromatic  waves  of  arbitrary 
wavelength.  The  results  are  compared  to  those  obtained  from 
Fourier  transformation  and  to  standard  deviations,  which 
highlight  the  characteristics  of  the  various  approaches  and 
the  advantages  afforded  by  the  MEM/HA  analysis. 

[11]  In  section  4,  we  combine  CRISTA  data  and  analysis 
to  isolate  enhanced  GW  temperature  amplitudes  observed  in 
the  stratosphere  over  southern  South  America  in  November 
1994.  The  high  spatial  resolution  in  all  three  dimensions 
provided  by  CRISTA  [Offermann  et  al,  1999;  Riese  et  al, 
1999;  Grossmann  et  al,  2002]  allows  us  to  investigate  the 
three-dimensional  structure  of  the  GWs.  These  investiga¬ 
tions  are  supported  by  state-of-the-art  meteorological  non¬ 
linear  mesoscale  model  simulations  that  act  to  “hindcast” 
wave  activity  here,  as  well  as  targeted  GW  ray-tracing 
studies.  The  detailed  comparisons  underline  the  necessity 
to  estimate  the  sensitivity  of  the  limb-sounding  technique 
for  understanding  the  observed  GW  stractures. 

2.  GW  Sensitivity  Estimation  for  an  Infrared 
Limb  Sounder 

[12]  For  the  interpretation  of  GW  data  obtained  from 
satellite  instruments,  detailed  knowledge  of  the  instrument’s 
sensitivity  to  disturbances  of  different  horizontal  and  verti¬ 
cal  wavelengths  is  essential  [Alexander,  1998;  Preusse  et  al, 
2000].  Wu  and  Waters  [1997]  and  McLandress  et  al.  [2000] 
provided  such  an  analysis  for  saturated  radiances  measured 
by  the  Microwave  Limb  Sounder  (MLS)  on  the  UARS 
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local  horizontal 


Figure  1.  Measuring  geometry  of  the  limb  sounder.  Two 
rays  with  different  tangent  heights  are  shown.  These  are 
obtained  by  tilting  the  primary  mirror  of  the  instrument. 


satellite.  For  optically  thin  conditions  one  might  estimate 
an  instrument’s  sensitivity  to  different  vertical  scales  of 
GWs  based  on  the  vertical  weighting  function  [Fetzer, 
1990].  Simply  speaking,  the  weighting  function  defines 
how  much  of  the  observed  radiation  stems  from  which 
part  of  the  atmosphere.  One  often  estimates  the  visibility  of 
a  small-scale  disturbance  by  convolving  the  disturbance 
field  with  the  weighting  function  (this  will  be  discussed 
later  in  this  section  in  more  detail).  To  estimate  the  sensi¬ 
tivity  to  different  horizontal  wavelengths  one  would 
analogously  invoke  the  horizontal  weighting  function, 
which  is  approximately  Gaussian.  This  gives  an  amplitude 
sensitivity 

=  (1) 

Ao 

where  Ag  is  the  amplitude  of  the  atmospheric  wave 
perturbation,  A  is  the  retrieved  amplitude,  is  the 
horizontal  wavelength  along  the  limb  direction  and  a  is 
the  width  of  the  horizontal  weighting  function  with  typical 
values  of  ^200  km.  Following  (1),  a  limb  sounder  has 
less  than  5%  sensitivity  to  a  500  km  wavelength  GW  and 
shorter  waves  would  vanish  completely  due  to  smearing 
along  the  ray  path.  Reasonable  sensitivity  is  found  only 
for  horizontal  wavelengths  longer  than  ^800  km.  This  low 
sensitivity  is  in  obvious  contradiction  to  the  mesoscale 
horizontal  structures  seen  in  temperature  data  from  the 
LIMS  and  the  CRISTA  instruments  [Fetzer,  1990;  Fetzer 
and  Gille,  1994;  Eckermann  and  Preusse,  1999;  Preusse 
et  al,  2001b]  indicating  that  (1)  does  not  adequately 
describe  the  sensitivity  of  these  instruments  to  GWs.  Thus, 
here  we  consider  a  two-dimensional  geometry  and 
generalize  the  results  to  three  dimensions  at  the  end  of 
the  section. 

2.1.  Limb-Sounding  Technique 

[i3]  The  limb-sounding  technique  has  been  described  in 
detail  elsewhere  [e.g.,  Bailey  and  Gille,  1986;  Gordley  et  al., 
1994;  Marshall  et  al.,  1994;  Riese,  1994;  Riese  et  al.,  1999] 
and  is  depicted  schematically  in  Figure  1.  The  satellite 


views  from  the  orbit  position  downward  through  the  atmos¬ 
pheric  limb.  The  instrument  measures  the  limb  radiance, 
produced  by  emission  and  reabsorption  along  the  ray  path. 
The  mathematical  description  of  this  radiative  transfer  is 
given  for  example  by  Gordley  et  al.  [1994].  Vertical 
distributions  are  obtained  by  scanning  the  beam  up  and 
down  in  the  atmosphere. 

[14]  Figure  2  shows  the  viewing  geometry  of  a  limb¬ 
sounding  instrument  projected  onto  a  simulated  GW 
temperature  oscillation.  Since  we  use  a  Cartesian  x-z 
coordinate  system,  the  curvature  of  the  Earth’s  surface 
now  shows  up  in  form  of  a  curved  ray  path.  If  we  assume 
that  the  atmosphere  is  optically  thin,  we  can  neglect  self¬ 
absorption  in  the  atmosphere  and  in  this  case  the  intensity 
measured  by  the  satellite  is  the  sum  of  all  emissions  along 
the  ray  path.  Atmospheric  densities  increase  as  we  follow 
the  ray  path  downward.  For  emitters  with  constant  mixing 
ratio  profiles  (e.g.,  CO2),  density  increases  proportionally 
as  the  total  air  density.  Thus,  the  radiance  emitted  from  a 
small  section  along  the  ray  is  proportional  to  the  atmos¬ 
pheric  density  multiplied  by  the  blackbody  radiance  for  the 
corresponding  atmospheric  temperature.  This  emitted  radi¬ 
ance  is  depicted  in  the  upper  panel  of  Figure  2  by  the 
Gaussian  curve  and  shows  the  weighting  that  has  to  be 
applied  to  an  atmospheric  temperature  fluctuation  in  order 
to  obtain  the  temperature  fluctuation  detected  by  the 
satellite  (cf.  (1)). 

[15]  The  two  ray  paths  shown  in  the  left  panel  of  Figure  2 
reach  their  lowest  altitude  at  the  tangent  point  at  x  =  500  km, 
where  the  weighting  function  peaks.  Also  given  in  the  upper 
panel  is  a  horizontal  projection  of  the  wave  oscillation, 
showing  the  horizontal  wavelength.  As  discussed  above,  a 
convolution  of  this  wave  structure  with  the  weighting 
function  would  lead  to  a  nearly  complete  smearing  out  of 
the  wave.  However,  we  can  observe  that  the  ray  path  is 
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Figure  2.  Viewing  geometry  of  a  limb  sounding  projected 
through  a  simulated  GW.  Right  panel  shows  the  temperature 
profile  used  as  model  input  (dashed)  and  the  modeled 
profile  at  the  tangent  point  (x  =  500  km;  solid  line).  The 
two-dimensional  temperature  fluctuations  and  two  limb  ray 
paths  are  shown  in  the  left  panel.  Top  panel  gives  the 
weighting  function  along  a  ray  path  and  the  horizontal 
waveshape  in  temperature  for  the  bottom  altitude.  For 
details,  see  text. 
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partly  coaligned  along  the  sloping  phase  front  of  the  GW  at 
X  -  200-400  km.  In  addition  the  phase  does  not  change 
greatly  between  x  -  400  km  and  the  tangent  point.  Thus,  to 
the  left  of  the  tangent  point  all  ray  sections  with  large  weight 
have  about  the  same  phase.  On  the  right,  phases  change 
rapidly  and  no  net  signal  can  be  expected  due  to  averaging 
along  this  limb  segment.  Thus,  the  phase  seen  between  x  - 
200  km  and  x  =  500  km  will  dominate  the  radiance  attributed 
to  this  tangent  height.  Considering  other  tangent  heights, 
different  phases  will  be  emphasized  and  a  vertical  wave 
oscillation  is  measured  by  the  satellite  instrument. 

2.2.  Analytical  Solution  of  the  Radiative  Transfer 

[i6]  We  can  quantify  all  this  by  solving  the  radiative 
transfer  equations  analytically  given  relevant  simplifica¬ 
tions.  The  analytical  solution  not  only  provides  a  better 
understanding  of  the  radiative  transfer  and  the  sensitivity  to 
GWs,  but  also  avoids  numerical  errors  due  to  grid  spacing. 
Thus  it  can  be  used  for  cross  checking  the  results  from  a  full 
numerical  model  of  the  measurement  and  retrieval  as  will  be 
done  in  section  2.3.  The  analytical  solution  we  will  develop 
uses  the  following  assumptions; 

1 .  The  blackbody  radiance  B(T)  can  be  linearly  expanded 
in  the  temperature  T  using  a  first  order  Taylor  series. 

2.  The  wavelength  region  considered  is  optically  thin. 

3.  The  background  atmosphere  is  isothermal. 

4.  The  wave  amplitude  remains  constant  with  altitude. 

[iv]  We  first  consider  the  radiance  emitted  along  the  ray 

path.  For  an  optically  thin  atmosphere  we  can  neglect  self¬ 
absorption  and  the  limb  radiance  becomes 

POO 

/=  /  t(s)B(T(s))ds,  (2) 


where  s  is  the  coordinate  along  the  ray  path.  The  emissivity 
e  along  a  small  part  of  the  ray  is  directly  proportional  to 
the  total  atmospheric  density  p.  Since  the  density  in  an 
isothermal  atmosphere  is  given  by  p  =  poe“^^^  we  can 
rewrite  this  equation  as 

/OO 

p,e-^I^B{T{s))ds.  (3) 

-OO 

Thus  W=  is  the  weight  which  has  to  be  applied  to  any 
atmospheric  temperature  perturbation  to  be  considered.  We 
now  express  this  weighting  function  in  terms  of  a  horizontal 
weighting  function  by  replacing  z  by  s\  Figure  3  shows  the 
geometry  of  a  single  tangent  ray.  In  contrast  to  Figure  2,  here 
the  coordinate  5  is  a  straight  line  and  the  coordinate  x  follows 
the  Earth’s  curvature.  From  Figure  3  we  obtain  the  relations  s' 
=  /i/  sin  (f)  —  =  V  and  ht  —  hi  costf)  ~  where 

hi  =  +  Zi  is  the  distance  from  the  Earth’s  center  to  the 

considered  ray  element  and  ht  =  Re  +  z,  is  the  distance  to  the 
tangent  point.  Since  z  =  zi  —  Zt'^Q  get: 


(4) 


Our  weight  takes  the  shape  of  a  Gaussian 


(5) 


Earth  center 


Figure  3.  Geometry  of  a  single  ray  path.  With  varying 
distance  s  from  the  tangent  point  the  altitude  z/  of  the  ray 
above  the  Earth’s  surface  and  the  distance  to  the  Earth’s 
center  hi  =  Rsarth  +  z/  changes.  Tangent  height  z^,  s,  and  z/  can 
be  all  connected  by  trigonometric  functions  in  tf)  (see  text). 


and  we  obtain  a  =  ^/HRe-  Choosing  H  =  6.5  km  and  Re  = 
6350  km  we  get  a  ~  200  km,  which  we  have  used  to 
evaluate  (1). 

[is]  We  now  calculate  the  deviation  of  the  intensity 
emitted  from  a  sinusoidally  perturbed  atmosphere  from 
the  intensity  for  the  background  atmosphere  by  applying 
the  weighting  function.  In  this  step  we  also  linearly  expand 
the  blackbody  radiation  as  a  function  of  temperature,  which 
requires  small  deviations  from  the  mean  temperature 

/OO 

fV(s)  A  sin(fa  -f  mz  +  'if}(z/))ds,  (6) 

-OO 

with  A  being  the  temperature  amplitude  and  k  =  and 
m  =  2t^/\  being  the  horizontal  and  vertical  wave  numbers, 
respectively.  Thus,  after  making  the  approximation  x  ~  5  the 
only  but  important  difference  between  (6)  and  the  assump¬ 
tions  leading  to  (1)  is  that  (6)  contains  a  height  dependence 
mz(s).  This  height  dependence  z(^)  is  defined  by  the 
constraint  that  the  air  volumes  considered  lie  along  the  limb 
ray.  The  dependence  of  the  wave  phase  on  tangent  height  is 
expressed  by  the  phase  r|)(Zi)  =  mz,  -I-  i(io-  For  simplicity  we 
assume  here  =  0  and  use  (4)  to  replace  z  by  ^  and  use 
addition  theorems  to  extract  from  the  integral.  In  addition, 
we  omit  integrals  of  antisymmetric  functions 

A/(z,)  oc 

A  cos(mzt)  cos(G)sin  (^m  exp  (^-  ds 

+A  sin(mz,)/!|^cos(G)cos(^m^)exp(^-24jyi  (7) 

An  analytical  solution  for  these  integrals  exists  \Gradshteyn 
and  Ryzhik,  1994].  The  relative  radiance  sensitivity  is 
obtained  by  dividing  the  intensity  variation  AI{Zf)  by  the 
mean  intensity  /(z,)  and  the  atmospheric  amplitude  A.  We 
therefore  normalize  to  the  weight 


W{s)  = 


(8) 
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Figure  4.  Estimated  sensitivity  of  an  infrared  limb 
sounder  as  obtained  using  an  analytical  approach.  Panel 
(a)  shows  the  radiance  amplitudes  as  percentage  deviations 
from  the  total  radiance  corresponding  to  an  atmospheric 
oscillation  with  a  1  K  temperature  amplitude.  Panel  (b) 
gives  the  ratio  of  retrieved  and  atmospheric  temperature 
amplitude  for  an  “idealized”  retrieval  without  any  restric¬ 
tions  from  the  instrument.  Panel  (c)  takes  sampling  effects 
into  account.  The  same  ratio  as  given  in  Panel  (b)  is  shown, 
but  as  expected  for  a  3  km  step-width  onion-peeling 
retrieval.  The  plots  are  based  on  analytical  expressions 
derived  in  the  text.  The  contour  spacing  is  0.1  in  Panel  (a) 
and  0.05  in  Panels  (b)  and  (c). 


and  take  into  account  a  factor  f  ||  stemming  from  the  linear 
expansion  of  the  blackbody  radiation.  We  finally  obtain  the 
sensitivity  of  the  total  ray  radiance  to  a  sinusoidal 
perturbation  in  the  atmosphere: 


SA,2D{H,k,m) 


1  dB  aI  c2 


cb^  \ 
4(c^  +a^))’ 


(9) 


where  a  =  m/2R^=  b  =  and  c  =  \I{2HR^. 

In  order  to  simulate  the  CRISTA  temperature  retrieval  we 
evaluate  at  12.6  p,/M  wavelength  of  the  IR  radiation  and  at 

a  mean  atmospheric  temperature  of  230  K.  The  phase  shift 
between  the  atmospheric  wave  in  temperature  and  the  wave 
in  the  limb  radiances  is 

.  1  c^b^ 

All)  =  - arctan  5 (10) 

^  2  \cJ  4(c2  +  a2)  ^  ’ 

[19]  Sa,2d{H,  k,  m)  is  plotted  in  Figure  4a  for  a  constant 
value  of  H  =  6.5  km.  We  find  that  in  general  short  vertical 
wavelengths  are  suppressed  by  the  radiative  transfer.  How¬ 
ever,  an  “idealized”  retrieval  would  recover  all  these 
degradations  in  the  vertical,  given  an  infinite  signal-to-noise 
ratio  of  the  instrument.  We  first  want  to  approximate  an 
idealized  retrieval,  which  does  not  further  degrade  the 
measured  signals  in  the  vertical,  because  this  allows  us  to 
separate  the  different  constraints  more  clearly.  Furthermore, 
it  allows  the  retrieval  step  width  to  be  changed  easily  (see 
below).  Some  retrieval  algorithms  take  into  account  hori¬ 
zontal  gradients  of  the  background  atmosphere.  However, 
the  fine-scale  horizontal  variations  of  a  mesoscale  GW 
cannot  be  taken  into  account  in  the  retrieval  process 
particularly  when  the  viewing  direction  is  not  aligned  with 
the  flight  vector.  Thus,  the  retrieval  will  recover  the  wave 
structure  as  if  it  were  measured  in  the  absence  of  any 
horizontal  structure,  i.e.,  for  infinite  horizontal  wavelengths. 
We  can  therefore  approximate  the  retrieval  by  dividing 
Sa,2d{H,  k,  m)  by  Sa.2d{H,  k  =  0,  m). 

[20]  Applying  this  “idealized”  retrieval  we  obtain 
Figure  4b.  From  this  figure  we  can  already  see  some  im¬ 
portant  features.  For  vertical  wavelengths  around  10  km  the 
limb  sounder  has  reasonable  sensitivity  (0.35),  even  for 
wavelengths  around  200  km.  This  is  much  better  than  the 
ID  estimations  from  (1).  We  also  see  that  for  short  horizon¬ 
tal  wavelengths  (<1000  km)  the  sensitivity  at  long  vertical 
wavelengths  is  diminished.  This  is  an  important  point  for  the 
intercomparison  of  satellite  climatologies  \Preusse  et  al, 
2000].  It  also  explains  the  low  sensitivity  predicted  by  (1). 
Setting  fl  =  0  in  (9)  we  get  (1).  Thus,  the  one-dimensional 
consideration  is  given  by  Sa.2d{H,  k,  m  =  0)  in  (9).  Figure  2 
for  this  case  would  show  vertically  aligned  phase  fronts  and 
the  limb  ray  can  never  be  oriented  parallel  to  the  phase 
fronts.  Thus  the  one-dimensional  estimation  (1)  is  a  worst 
case,  which  is  never  found  in  reality. 

[21]  From  Figure  4b  we  would  expect  to  see  GWs  with 
short  vertical  wavelengths  best.  However,  this  was  obtained 
assuming  an  “idealized  retrieval”,  whereas  real  retrievals 
have  additional  limitations.  The  smallest  resolvable  vertical 
scale  visible  in  the  CRISTA  data  is  limited  by  the  vertical 
sampling  and  the  retrieval  process.  During  the  first  flight  in 
November  1994,  CRISTA  sampled  with  a  vertical  resolu¬ 
tion  of  1.5  km  [Riese  et  al,  1999].  However,  except  for 
some  special  observational  modes,  we  have  divided  the  data 
set  into  two  subsets  of  3.0  km  vertical  step  width,  which  are 
retrieved  separately  and  reemerged  afterward  {Riese  et  al, 
1999].  Since  a  simple  onion-peeling  scheme  is  used  for  the 
retrieval  {Riese,  1994;  Riese  et  al,  1999],  temperature  and 
trace-gas  mixing  ratios  are  interpolated  linearly  between  the 
altitudes  considered.  Thus,  the  retrieval  linearly  interpolates 
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throughout  any  altitude  structure  which  is  smaller  than  the 
retrieval  step  width  of  3.0  km. 

[22]  We  simulate  this  in  the  analytical  approach  by 
applying  a  boxcar  filter  of  the  size  of  the  retrieval  step 
width  Az.  The  smoothing  inherent  in  the  retrieval  process  is 
then  given  by 


^(X,) 


211  Az 


(11) 


Applying  this  “cut  off”  due  to  the  retrieval  step  to  our 
“idealized  retrieved”  sensitivity  field  yields  Figure  4c.  This 
results  in  a  realistic  decrease  in  sensitivity  at  short  vertical 
wavelengths. 

2.3.  High  Spatial  Resolution  Radiative  Transfer 
Calculations 

[23]  The  analytical  approach  is  based  on  several  approx¬ 
imations  which  may  limit  the  general  applicability  of  the 
results.  A  more  general  approach  is  provided  here  by  a 
numerical  modeling  that  is  divided  into  three  steps.  In  the 
first  step  a  representative  background  atmosphere  is  per¬ 
turbed  using  a  linear  two-dimensional  (2D)  GW  oscillation 
based  on  WKB  theory.  Next  radiative  transfer  simulations 
throughout  the  model  atmosphere  are  performed  leading  to 
vertical  radiance  profiles  such  as  those  measured  by 
CRISTA.  (The  simulations  are  idealized  in  a  sense  that 
instrumental  error  sources  such  as  radiance  errors  or  noise 
are  not  taken  into  account.)  The  simulated  input  amplitudes 
are  chosen  to  be  typically  a  few  Kelvin  and  will  be  easily 
identified  by  the  MEM/HA  given  a  reasonable  sensitivity 
value.  The  third  step  involves  retrieval  and  data  analysis, 
which  are  performed  following  the  CRISTA  standard  data 
evaluation. 

[24]  The  model  used  is  a  two-dimensional,  linear  WKB 
model  for  MWs  [Schoeberl,  1985].  Horizontal  wavelength 
and  altitude  profiles  of  wind  speed  and  potential  temper¬ 
ature  are  input  to  the  model.  For  our  calculations  we  have 
chosen  a  wind  profile  that  is  constant  with  altitude.  The 
value  of  the  horizontal  wind  speed  U  sets  the  vertical  wave 
number  m 


where  N  is  the  buoyancy  frequency,  k  is  the  horizontal  wave 
number  and  H  is  the  scale  height.  The  bottom  layer  of 
potential  temperature  follows  the  surface  topography,  which 
here  is  sinusoidal.  For  all  other  altitudes  the  displacement  of 
the  isentropes  is  calculated  according  to  linear  wave  theory. 
Amplitude  growth  with  altitude  follows  the  WKB  relation 
[e.g.,  Lindzen,  1981] 


data  extrapolated  to  the  ground  and  to  about  80  km.  The 
extrapolation  to  80  km  was  necessary  for  the  forward 
radiance  calculations. 

[25]  The  primary  model  output  is  a  perturbed  potential 
temperature  field,  whereas  the  radiative  transfer  calculations 
are  based  on  temperatures,  pressures  and  trace  gases  on  a 
geometric  altitude  grid. 

[26]  Therefore  the  atmosphere  is  rebuilt  hydrostatically 
from  the  perturbed  potential  temperature  field  and  temper¬ 
ature  and  pressure  are  evaluated.  The  results  are  stored  on  a 
grid  of  10  km  in  the  horizontal  and  500  m  in  the  vertical.  It 
should  be  noted  that  ozone,  which  is  the  most  important 
interfering  trace  species  in  the  retrieval  process,  is  perturbed 
by  the  wave  also. 

[27]  From  this  model  atmosphere  we  obtain  altitude 
profiles  of  limb  radiances  as  would  be  measured  by 
CRISTA  with  the  use  of  radiative  transfer  calculations  based 
on  the  BANDPAK  radiance  model  [Marshall  et  al,  1994; 
Gordley  et  al.,  1994].  Every  combination  of  horizontal  and 
vertical  wavelength  had  to  be  simulated  individually.  For 
each  profile  we  define  a  vertical  grid  of  tangent  heights  with 
a  spacing  of  200  m,  oversampling  the  1.5  km  vertical 
sampling  of  the  CRISTA  instrument.  This  allows  us  to 
convolve  the  radiances  resulting  from  an  infinitesimally 
thin  ray  path  with  the  vertical  field  of  view  of  the  CRISTA 
instrument.  For  each  ray  path  calculation,  the  atmosphere  is 
divided  into  100  layers.  These  layers  define  the  spatial 
resolution  of  the  radiative  transfer  calculations.  The  vertical 
spacing  of  the  layers  is  dense  immediately  above  the  tangent 
height  and  increases  with  altitude.  This  means  highest 
spatial  resolution  is  obtained  where  the  weighting  function 
peaks.  A  fine  vertical  resolution  of  the  layers  is  especially 
important  for  short  horizontal  wavelengths  Xj^,  because  in 
this  case  the  phase  of  the  wave  changes  rapidly  in  the  x- 
direction,  but  close  to  the  tangent  point  the  height  of  the  ray 
path  varies  slowly  with  x  (see  Figure  2).  It  turned  out  that  in 
some  cases  a  vertical  resolution  of  around  30  m  is  necessary 
for  the  layers  immediately  above  the  tangent  point.  At  the 
altitude  of  each  layer  the  exact  position  of  the  ray  path  is 
detennined  and  the  simulated  atmospheric  parameters  given 
by  the  2D  model  are  interpolated  to  these  ray  positions.  For 
these  radiative  transfer  calculations  none  of  the  approxima¬ 
tions  necessary  for  the  analytical  approach  are  made. 

[28]  The  simulated  radiance  profiles  are  retrieved  using 
the  standard  CRISTA  retrieval  algorithms,  including  all 
smoothing  procedures  [Riese  et  al,  1999].  Finally  the 
retrieved  temperature  profiles  are  analyzed  for  GWs  employ¬ 
ing  the  MEM/HA  method  described  in  section  3.  The 
resulting  amplitudes  and  wavelengths  can  be  compared  to 
the  input  values.  To  determine  the  amplitude  sensitivity  the 
output  of  the  simulated  retrieval  is  compared  to  profiles 
obtained  by  vertical  sections  throughout  the  GW  model 
output  at  the  tangent  point. 


w(z)  =  vv(0), (13) 
y  "((z) 

where  w  is  the  vertical  velocity  amplitude.  In  addition, 
calculations  were  perfonned  where  the  density  stratification 
term  e^^^^  was  omitted.  The  background  potential  tempera¬ 
ture  is  taken  from  a  zonal  mean  cross  section  of  CRISTA 


2.4.  Results  of  the  Numerical  Simulatious 

[29]  Visual  inspection  of  simulated  altitude  profiles 
reveals  a  phase  shift  between  the  wave  signatures  in  temper¬ 
atures  extracted  directly  from  the  2D-model  solution  and 
those  obtained  by  simulated  radiative  transfer  and  retrieval. 
These  phase  shifts  6([)(z)  depend  on  the  background  atmos¬ 
phere  and  hence  on  the  altitude.  This  can  be  noted  from 
(10),  where  Aij)  depends  on  the  scale  height  H  via  c. 
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Figure  5.  Scatterplot  of  retrieved  versus  theoretical  vertical  wavelength.  The  solid  line  indicates  a 
perfect  retrieval  of  the  vertical  wavelength.  The  observed  scatter  of  the  singular  points  as  retrieved  from 
simulated  CRISTA  profiles  is  due  to  radiative  transfer  effects  and  the  analysis  method.  The  scatter  gives  a 
good  indication  of  the  precision  of  the  vertical  wavelength  retrieved  from  real  CRISTA  measurements. 


To  demonstrate  the  way  in  which  height  dependent  phase 
shifts  can  alter  the  vertical  wavelength,  we  consider  a  sinus¬ 
oidal  wave  in  a  single  retrieved  profile  T'  =  T  sin(mz  +  (|)  + 
64)(z))  =  T  sin  i[),  with  tf)  being  the  phase  at  z  =  0.  For 
example,  if  the  phase  shift  6(|)(z)  in  a  certain  altitude  range 
decreases,  the  total  phase  of  the  wave  changes  more 
slowly  than  we  would  expect  from  the  vertical  wave  number 
m.  This  results  in  an  apparently  longer  vertical  wavelength 
in  the  retrieved  temperatures  and  by  the  same  token  an 
increasing  phase  shift  results  in  a  shorter  wavelength.  This 
variation  of  the  vertical  wavelength  is  shown  in  Figure  5. 
The  vertical  wavelength  deduced  from  the  retrieved  temper¬ 
atures  is  compared  to  the  input  wavelength  from  (12).  We 
find  a  scatter  of  ^20%  in  retrieved  wavelengths. 

[3o]  The  simulations  employed  to  produce  Figure  5  are 
based  on  a  zonal  mean  background  atmosphere  and  there¬ 
fore  embody  realistic  geographical  variations.  In  addition, 
the  same  method  as  used  here  is  also  used  to  determine  the 
vertical  wavelength  of  GWs  observed  in  real  flight  profiles 
(see  section  3).  Thus,  the  scatter  observed  in  Figure  5 
represents  physical  effects  from  the  radiative  transfer  as 
well  as  influences  of  the  analysis  method  and  is  a  good 


estimate  of  the  precision  of  the  vertical  wavelength  retrieved 
from  CRISTA  temperatures. 

[3i]  GW  simulations  were  performed  for  constant  ampli¬ 
tudes  and  amplitudes  growing  exponentially  with  height.  In 
the  real  atmosphere  both  cases  are  likely  to  be  observed. 
Exponential  growth  is  expected  where  GWs  propagate 
conservatively  without  interaction  with  the  background 
flow.  A  constant  amplitude  reflects  the  classical  case 
described  by  Lindzen  [1981],  where  GWs  are  breaking  and 
retain  their  saturation  amplitude.  It  should  further  be  noted 
that  even  in  the  case  of  exponential  growth  the  amplitude 
increases  with  twice  the  scale  height  of  the  density  decrease. 
Thus,  the  weight  along  the  ray  path  decreases  faster  than  the 
wave  amplitude  increases  and  the  total  contribution  to  the 
measured  radiance  decreases  exponentially  with  altitude. 
Results  from  growing  amplitude  as  well  as  constant  ampli¬ 
tude  simulations  are  shown  in  Figure  6.  Panel  (a)  reproduces 
Figure  4c.  Panels  (b)  and  (c)  show  results  for  constant 
amplitude  and  exponentially  increasing  amplitude  waves, 
respectively.  The  values  of  the  numerical  simulations  are  an 
average  over  altitudes  of  25-60  km  and  13  latitude  bins. 
This  average  increases  the  statistical  significance  of  the 
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Figure  6.  Estimated  sensitivity  of  retrieved  CRISTA 
temperatures  to  GW  oscillations.  Radiative  transfer  as  well 
as  retrieval  effects  are  included.  Panel  (a)  reproduces  Figure 
4c.  Panels  (b)  and  (c)  show  the  results  of  the  numerical 
simulations  with  constant  GW  amplitudes  and  amplitudes 
growing  exponentially  with  altitude,  respectively.  All  three 
figures  exhibit  the  same  salient  features,  particularly  the 
strong  decrease  of  sensitivity  at  around  200  km  horizontal 
wavelength  and  the  decreased  sensitivity  to  long  vertical 
wavelength  for  a  given  horizontal  wavelength.  Panels  (b) 
and  (c)  are  very  similar  in  the  horizontal  wavelength  range 
from  300  to  800  km,  allowing  us  to  deduce  a  one¬ 
dimensional  correction  factor,  which  can  be  applied  to 
measured  CRISTA  data. 


given  values.  The  three  panels  have  the  same  salient  features, 
specifically  a  strong  decrease  of  sensitivity  at  \  ^  200  km 
and  a  decrease  of  sensitivity  at  long  vertical  wavelengths  for 
a  given  horizontal  wavelength.  The  results  of  the  full 


numerical  modeling  (Panels  (b)  and  (c))  reveal  an  “island” 
of  high  sensitivity  at  \  «  350  km  and  «  6-8  km  vertical 
wavelength,  which  is  not  found  in  the  analytical  solution. 
For  horizontal  wavelengths  between  300  and  800  km. 
Panels  (b)  and  (c)  agree  reasonably  and  exhibit  only  weak 
dependence  on  the  horizontal  wavelength.  Thus,  a  one¬ 
dimensional  sensitivity  function  averaged  over  the  300- 
800  km  horizontal  wavelength  range  and  depending  solely 
on  the  vertical  wavelength  is  representative  of  the  waves 
contained  in  the  residual  temperatures.  This  is  an  impor¬ 
tant  point,  because  it  allows  us  to  correct  amplitudes 
retrieved  from  CRISTA  measurements  for  visibility  effects 
noted  in  Figures  4  and  6  [see  also  Preusse  et  al,  2000]. 
This  would  not  be  possible  if  a  strong  dependence  on 
horizontal  wavelength  was  present,  since  it  is  often  diffi¬ 
cult  to  deduce  the  horizontal  wavelengths  as  will  be  seen 
in  section  4.  We  will  discuss  in  section  3.3  to  what  extent 
the  assumption  is  valid  that  horizontal  wavelengths  are 
shorter  than  1000  km. 

[32]  It  should  be  noted  that  there  is  an  artificial  back¬ 
ground  of  relatively  high  sensitivity  (0.2)  at  very  short 
horizontal  wavelengths.  This  is  spurious  and  due  to  a 
vertical  spacing  of  the  forward  calculations  that  is  still 
too  coarse.  Inspection  of  single  profiles  indicated  that  the 
altitude  layering  has  to  be  even  finer  than  the  30  m 
spacing  we  have  used  for  these  simulations.  Here  the 
analytical  approach,  which  is  not  restricted  by  numerical 
errors  and  grid  spacing,  supports  the  finding  of  vanishing 
sensitivity  to  the  very  short  horizontal  scales  less  than 
about  200  km.  For  longer  horizontal  wavelengths  the 
numerical  modeling  can  be  expected  to  yield  reliable 
results. 

[33]  So  far  we  have  discussed  sensitivity  variations 
versus  wavelength.  However,  an  important  question  is 
whether  the  observational  technique  itself  can  cause  artifi¬ 
cial  geographical  structures  in  the  observed  GW  variances. 
Thus,  we  now  keep  the  horizontal  and  vertical  wavelength 
fixed  by  averaging  all  values  in  a  box  of  X^,  =  300-800  km 
and  Xz  =  8-12  km.  Figure  7a  plots  these  values  versus 
latitude  and  altitude.  A  second  wavelength  box  of  X^:  = 
300-800  km  and  Xz  =  15-20  km  results  in  the  distribution 
shown  in  Figure  7b.  It  should  be  noted  that  the  MEM/HA 
method  results  in  wave  amplitudes  which  are  independent 
of  the  phase  of  the  wave  relative  to  the  analysis  interval.  In 
addition,  the  average  of  X^,  =  300-800  km  and  different 
vertical  wavelengths  is  also  an  average  over  different 
phases  of  the  waves.  The  observed  stmctures  are  therefore 
reliable  and  not  due  to  a  specific  phase  choice  made  when 
initializing  the  simulations.  The  observed  geographical 
dependence  is  small  for  the  short  vertical  wavelengths 
(Panel  a)  and  varies  by  up  to  ±30%  for  the  longer  vertical 
wavelengths  shown  in  Panel  (b).  However,  this  higher 
value  corresponds  to  about  2  dB  for  squared  amplitudes 
and  we  will  show  in  section  3  and  the  companion  paper 
(Preusse  et  al.,  in  preparation,  2002)  that  these  variations 
are  small  compared  to  the  observed  geophysical  variations 
of  GW  variances. 

[34]  One  last  point  is  whether  there  are  changes  in 
the  phase  of  the  observed  waves.  The  analytical  ap¬ 
proach  indicates  that  the  radiative  transfer  (10)  as  well  as 
the  retrieval  induce  phase  shifts.  These  are  expected 
mostly  where  the  radiance  sensitivity  is  small,  namely  for 
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Figure  7.  Dependence  of  the  sensitivity  on  the  location  of 
the  measurement.  Short  vertical  wavelengths  (V  =  8- 
12  km)  given  in  Panel  (a)  depend  only  weakly  on  altitude 
and  are  nearly  independent  of  latitude.  The  geographical 
dependence  of  the  longer  vertical  wavelengths  (Xz  =  15- 
20  km)  is  stronger.  However,  as  shown  in  sections  3  and  5, 
the  observed  variations  in  GW  activity  are  much  stronger. 

short  vertical  and  short  horizontal  wavelengths  (compare 
Figures  4  and  6).  Results  from  numerical  simulations  are 
shown  in  Figure  8,  which  plots  the  absolute  values  of  the 
phase  differences  between  2D-model  wave  input  and 
retrieved  temperature  profiles.  These  phase  shifts  are  small 
for  most  combinations  of  horizontal  and  vertical  wave¬ 
length.  The  phase  shifts  therefore  do  not  prevent  us  from 
deriving  horizontal  wave  structures  from  the  CRISTA 
temperature  data. 

2.5.  Generalization  to  Three  Dimensions 

[35]  The  sensitivity  estimates  shown  in  Figures  4  and  6 
are  based  on  two-dimensional  radiative  transfer  calcula¬ 
tions.  However,  the  real  atmosphere  is  three  dimensional 
and  it  is  likely  that  the  instrument  views  waves  at  various 
angles.  This  is  illustrated  in  Figure  9,  which  shows  the 
horizontal  projection  of  a  limb  ray  and  the  horizontal  phase 
fronts  of  a  GW.  The  horizontal  wavelength  along  the  line  of 
sight  (LOS)  is  a  factor  l/(cos  a)  larger  than  the  real 
horizontal  wavelength  X^..  Thus,  for  waves  viewed  at  an 
angle  a  to  the  phase  fronts  we  have  to  scale  the  x-axes  of 
Figures  4  and  6  by  cos(a).  Taking  a  =  60°  as  an  example, 
the  shortest  horizontal  wavelength  X^  for  which  the  sensi¬ 
tivity  is  larger  than  0.5  is  now  100-120  km  instead  of  200- 
240  km  in  the  two-dimensional  calculations  (Figures  6b 
and  6c),  which  tacitly  assume  a  =  0. 


Figure  8.  Phase  shift  due  to  the  radiative  transfer.  The 
values  are  mostly  small  except  for  wavelengths,  for  which 
the  small  sensitivities  have  been  found  in  Figure  6. 


[36]  We  can  assume  that  the  case  where  the  LOS  is 
coaligned  with  the  horizontal  wave  vector  is  the  exception 
and  that  nonzero  view  angles  a  frequently  occur.  Thus,  a 
fraction  of  the  atmospheric  waves  with  horizontal  wave¬ 
lengths  between  100  km  and  200  km  is  observed  by  the 
limb  viewer.  How  large  this  fraction  is  may  depend  on  a 
possible  general  preference  of  the  wave  propagation  direc¬ 
tion  and  the  actual  view  direction  of  the  instrument.  This  is 
discussed  in  some  detail  in  section  4.  However,  for  very 
large  view  angles  a  this  relation  may  not  hold  any  longer 
because  this  would  assume  infinite  plane  horizontal  wave 
fronts. 

[37]  Employing  the  one-dimensional  correction  function 
[Preusse  et  al,  2000]  one  can  rescale  the  measured  GW 
amplitudes  to  estimate  the  “true”  atmospheric  amplitudes 
of  mesoscale  GWs  with  vertical  wavelengths  \  ^  5  km. 
The  decrease  of  sensitivity  at  the  short  horizontal  wave- 


Figure  9.  Horizontal  projection  of  an  arbitrarily  oriented 
GW  and  the  line  of  sight  (LOS).  Instead  of  viewing 
perpendicular  to  the  wave  fronts  (gray  shading),  there  is  an 
angle  a  between  the  line  of  sight  and  the  horizontal  wave 
vector  k. 
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length  edge  is  very  steep  and  we  therefore  can  assume  that  a 
wave  is  either  scaled  correctly  or  not  measured  at  all.  The 
importance  of  rescaling  the  measured  amplitudes  is  pointed 
out  by  Preusse  et  al.  [2000], 

3.  Data  Acquisition  and  Analysis 

[38]  We  have  shown  in  the  previous  section  that  meso- 
scale  GWs  cause  detectable  signatures  in  temperature  alti¬ 
tude  profdes  acquired  by  an  infrared  emission  limb  sounder. 
Such  vertical  temperature  profiles  were  sampled  by  the 
CRISTA  instrument  with  high  spatial  resolution  and  nearly 
global  coverage.  In  this  section  we  briefly  introduce  the 
instrument  and  discuss  the  mathematical  tools  used  to 
analyze  the  data  to  extract  the  vertical  wavelength,  ampli¬ 
tude  and  phase  of  GW  oscillations  in  acquired  temperature 
profiles. 

3.1.  Measurement  Technique  and  Detrending 

[39]  The  CRISTA  instrument  was  designed  to  acquire 
very  high  spatial  resolution  data.  For  this  purpose  CRISTA 
uses  three  telescopes  viewing  18°  apart  to  measure  height 
profiles  of  thermal  emission  limb  spectra  in  the  middle  and 
far  infrared  at  a  spectral  resolution  of  about  500  [Offermann 
et  al,  1999].  The  instrument  was  successfully  flown  during 
two  space  shuttle  missions,  STS  66  (November  1994)  and 
STS  85  (August  1997),  here  referred  to  as  CRISTA- 1  and 
CRISTA-2,  respectively.  Detailed  information  on  the  first 
mission  can  be  found  in  the  works  of  Offermann  et  al 
[1999]  and  Riese  et  al  [1999].  The  second  mission  is 
described  by  Grossmann  et  al  [2002].  On  both  flights  the 
spacecraft  orbited  at  300  km  altitude  with  a  57°  inclination. 
This  results  in  ^16  orbits  a  day.  Observations  range  from 
^57°S  to  ^65°N  and  from  ^72°S  to  ^72°N  during 
CRISTA- 1  and  CRISTA-2,  respectively. 

[40]  The  retrievals  of  stratospheric  temperatures  from 
CO2  band  emissions  around  12.6  qm  and  mesospheric 
temperatures  from  the  15  qm  CO2  band  are  described  by 
Riese  et  al  [1999].  The  vertical  sampling  steps  are  1.5  km 
and  2.0  km  for  CRISTA-1  and  CRISTA-2,  respectively. 
Given  a  sampling  time  of  ^1  s  per  spectrum,  a  complete 
altitude  profile  consisting  of  20  to  43  altitude  steps  is 
obtained  in  25  to  50  s  depending  on  the  measuring  mode. 
The  lower  boundary  of  the  stratospheric  measurements  was 
15  km  for  CRISTA- 1.  The  upper  boundary  depends  on 
measuring  mode.  For  an  overview  of  the  different  measur¬ 
ing  modes,  see  the  works  of  Riese  et  al  [1999]  and 
Grossmann  et  al  [2002]  for  CRISTA- 1  and  CRISTA-2 
observing  modes,  respectively.  In  the  stratosphere  and 
lower  mesosphere  the  precision  of  the  retrieved  temperature 
data  was  estimated  to  be  better  than  0.4  K  for  both 
missions. 

[41]  To  obtain  residual  temperature  fluctuations,  ret¬ 
rieved  temperature  data  were  detrended  by  subtracting  a 
background  atmosphere  estimated  by  a  Kalman  filter. 
Common  modes  of  planetary  waves  in  the  stratosphere 
are  expected  to  have  zonal  wave  numbers  up  to  4  and 
inspection  of  CRISTA  data  shows  no  significant  contribu¬ 
tions  from  higher  wave  numbers  (Oberheide,  personal 
communication,  2001).  Thus  a  0-6  wave  number  Kalman 
filter  [Rodgers,  1977;  Ern,  1993]  should  account  for  all 
planetary  wave  structures  in  the  background  [Fetzer  and 


Figure  10.  Example  of  the  combined  MEM  and  HA.  The 
left  panel  shows  an  altitude  profile  of  measured  CRISTA 
temperature  residuals  (solid)  and  the  reconstruction  from  the 
fit  (dashed).  The  profile  was  acquired  on  9  November, 
eastward  of  Chile  at  24°  S  and  82°  W.  The  right  panel  plots 
its  detrended  MEM  spectrum,  the  spectral  peaks  (bottom  to 
top  vertical  lines),  and  the  two  leading  amplitudes  (solid  and 
hatched  bars).  Spectral  density  scale  is  logarithmic. 


Gille,  1994;  Preusse  et  al,  2001a].  In  addition,  ascending 
and  descending  orbit  nodes  are  treated  separately  to  take 
into  account  tidal  structures,  which  become  important  in 
the  mesosphere  [Ward  et  al,  1999;  Oberheide  et  al, 
2000]. 

[42]  The  Kalman  filter  is  estimated  on  a  regular  altitude 
grid  of  1.5  km  spacing.  The  coefficients  are  estimated 
for  latitude  bands  of  0.5°  and  smoothed  afterward  with  a 
7-point  running  mean  in  latitude.  The  application  of  the 
Kalman  filter  to  the  CRISTA  data  is  described  in  detail  by 
Ern  [1993].  A  comparison  of  nonfiltered  and  Kalman- 
filtered  CRISTA  data  was  presented  by  Offermann  et  al 
[1999]. 

[43]  The  background  temperature  estimated  by  the  Kal¬ 
man  filter  is  interpolated  in  altitude  and  latitude  to  the 
precise  measurement  locations  of  the  individual  height 
profiles.  It  should  be  noted  that  the  estimation  of  the 
Kalman  filter  coefficients  is  performed  on  each  altitude 
layer  separately.  Therefore  no  degradation  of  the  vertical 
resolution  of  local  fluctuations  should  result  from  this 
detrending. 

3.2.  Analysis  Method 

[44]  Height  profiles  of  Kalman  detrended  residual  temper¬ 
ature  fluctuations  are  analyzed  by  a  combination  of  MEM 
[Press  et  al,  1992]  and  HA  as  illustrated  in  Figure  10.  The 
left  panel  shows  the  Kalman  detrended  vertical  temperature 
profile  as  the  solid  curve.  This  profile  is  spectrally  analyzed 
using  MEM.  The  MEM  is  performed  on  the  complete  height 
profile,  but  can  also  identify  wave  contributions  visible  only 
at  subintervals  of  the  profile.  A  spectrum  obtained  by  MEM 
generally  has  a  red  noise  contribution  which  particularly 
interferes  with  the  spectral  components  of  interest,  if  the 
number  of  data  points  in  a  given  analyzed  data  set  is  low. 
Thus,  weak  or  broad  peaks  may  be  difficult  to  detect  in  the 
presence  of  the  red  noise  background.  Therefore  we  here 
detrend  the  MEM  spectrum  by  a  quadratic  fit  through  the 
minima  of  the  spectrum.  This  detrended  MEM  spectrum 
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is  given  in  the  right  panel  of  Figure  10.  Thin  dashed 
vertical  lines  indicate  the  maxima  in  the  spectrum  as 
identified  by  a  peak  finding  algorithm.  These  wavelength 
peaks  are  the  basis  of  the  HA  performed  on  a  sliding 
interval  shifted  in  sfeps  of  one  sampling  point  through  the 
profile.  In  an  iterative  process  the  points  in  the  sliding 
interval  are  fitted  by  sinusoids  and  the  wave  with  the 
largest  amplitude  is  removed  from  the  data.  The  two 
largest  amplitudes  so  obtained  and  their  corresponding 
wavelengths  and  phases  are  attributed  to  the  central 
altitude  of  the  sliding  interval.  It  should  be  noted  that 
in  the  case  of  double  peaks,  which  can  occur  if  an  MEM 
is  applied  to  a  small  data  set  {Press  et  al,  1992],  all  the 
wave  amplitude  is  attributed  to  one  of  the  two  peaks.  This 
is  due  to  the  subtraction  of  the  fitted  wave  from  the  data 
before  applying  the  next  fit.  Hence  the  second  largest 
amplitude  stored  for  further  investigations  will  contain  an 
independent  wave  contribution. 

[45]  The  solid  dark  and  hatched  bars  in  the  right  panel  of 
Figure  10  indicate  the  wavelength  belonging  to  the  strongest 
and  second  strongest  amplitude  in  the  profde  on  the  left, 
respectively.  For  example,  the  two  temperature  maxima  at 
around  25  km  and  35  km  altitude  correspond  to  a  dominat¬ 
ing  wavelength  of  about  10  km  in  the  lower  altitudes.  The 
major  contribution  at  higher  altitudes  has  a  very  long 
wavelength  {k,2\  km)  visible  as  a  very  wide  trough  at 
around  45  km.  Above  43  km  the  second  wave  component 
captures  a  short  wavelength  contribution  of  around  5  km 
vertical  wavelength.  The  dashed  curve  in  the  left  panel 
shows  the  profile  reconsfrucfed  from  the  two  wave  compo¬ 
nents  with  the  strongest  amplitudes,  evaluated  at  each 
altitude  step  from  the  MEM/HA  values.  It  can  be  seen  that 
the  profde  is  locally  well  described  by  two  wave  compo¬ 
nents  which  may  vary  with  altitude. 

[46]  The  example  given  in  Figure  10  is,  however,  not 
really  typical.  Normally  the  altitude  profiles  look  less  rich 
and  only  fwo  wavelengths  are  identified  by  the  MEM. 
These  profiles  can  be  even  more  accurately  described  by 
the  MEM/HA  decomposition. 

[47]  Since  for  each  amplitude  value  the  corresponding 
vertical  wavelength  is  known,  the  results  can  be  corrected 
for  degradation  due  to  retrieval  and  radiative  transport 
following  the  method  described  in  section  2.4. 

3.3.  Method  Intercomparison 

[48]  The  MEM/HA  analysis  method  assumes  that  two 
monochromatic  waves  specify  almost  completely  the 
spectral  distribution  of  wave  energy  in  the  vertical  direc¬ 
tion  in  every  profile.  For  two-dimensional  MWs  this  is 
often  a  reasonable  assumption,  because  the  vertical  wave¬ 
length  is  governed  by  the  buoyancy  frequency  and  the 
wind  speed  only  [Queney,  1948;  Eckermann  and  Preusse, 
1999].  In  addition,  there  are  indications  from  lidar  meas¬ 
urements,  rocket  soundings,  radiosondes  and  radar  meas¬ 
urements  that  a  few  monochromatic  wave  components  are 
sufficient  to  describe  the  observed  temperature  fluctuations 
[e.g..  Cot  and  Barat,  1986;  Riggin  et  al.,  1995;  Sica  and 
Russell,  1999;  Bittner  et  al.,  1997].  However,  we  also 
investigated  the  general  applicability  of  the  method  for 
satellite  GW  analysis  on  a  statistical  basis.  Here  we 
compare  the  results  of  the  MEM/HA  to  those  obtained 
by  a  Fourier  Transfonn  (FT)  on  a  sliding  window  and  to 


Figure  11.  Leakage  functions  for  a  13-point  Fourier 
transform  at  1.5  km  vertical  sampling.  The  longest  vertical 
wavelength  is  about  20  km  and  the  shortest  wavelength  is 
20/6  «  3.3  km. 


conventional  standard  deviations,  the  latter  as  used  by 
Fetzer  and  Gille  [1994]. 

[49]  While  the  length  of  the  sliding  window  of  the 
MEM/HA  can  be  chosen  almost  arbitrarily,  the  length  of 
the  sliding  window  used  for  the  FT  is  directly  coupled  to 
the  wavelengths  resolved  by  the  FT.  For  this  comparison 
we  use  a  13-point  window  which  is  equivalent  to  an 
interval  length  of  ~20  km,  given  a  1.5  km  sampling  in 
altitude.  To  suppress  leakage,  a  Welch  window  is  applied 
to  the  data  before  the  FT.  The  spectral  power  is  nonnalized 
utilizing  Parseval’s  theorem  [Press  et  al.,  1992].  Figure  11 
shows  the  leakage  functions  of  the  different  Fourier 
components  dependent  on  the  vertical  wavelength  calcu¬ 
lated  according  to  Press  et  al.  [1992].  These  functions 
quantify  the  leakage  between  the  several  wavelength  bins, 
i.e.,  they  show  the  magnitude  of  the  contribution  from 
other  wavelengths  to  the  considered  wavelength  bin.  The 
wavelengths  of  maximum  sensitivity  (centers  of  the  wave¬ 
length  bins)  are  given  by  integer  fractions  of  the  length  of 
the  analysis  interval. 

[50]  We  first  compare  the  results  of  the  Fourier  transform 
to  the  zonal  mean  standard  deviations  calculated  from 
Kalman-detrended  temperature  fluctuations.  Since  Parse¬ 
val’s  theorem  has  been  used  to  normalize  the  results  of 
the  FT  for  a  stationary  process  we  would  expect  the  stand¬ 
ard  deviations  to  agree  with  the  sum  of  the  spectral  energy 
of  all  frequency  bins.  In  Figure  12  we  compare  the  standard 
deviations  (Panel  a)  with  the  sum  of  spectral  energy 
obtained  from  the  FT  method  (Panel  b).  Both  are  given  in 
dB  over  1  ranging  from  0  to  15  dB.  Panel  (c)  gives  the 
ratio  between  the  standard  deviation  and  the  FT  variances. 
At  altitudes  higher  than  about  30  km  the  patterns  in  Panels 
(a)  and  (b)  are  very  similar.  Thus,  on  average  the  ratio  is 
close  to  1  and  small-scale  patterns  range  from  0.7  to  1.2. 
Around  25  km  we  observe  that  the  FT  values  are  much 
higher.  This  is  due  to  the  vertical  width  of  the  sliding 
analysis  window.  A  value  obtained  by  the  FT  method  that 
is  attributed  to  25  km  contains  contributions  from  a  band  of 
altitudes  between  15  km  and  35  km.  In  particular  the  lower 
boundary  here  is  problematic  because  of  the  very  enhanced 
tropical  maximum  in  the  standard  deviations  observed 
between  15  km  and  20  km  altitude  (Figure  12a).  The  FT 
projects  this  maximum  to  higher  altitudes.  In  addition,  the 
subtropical  minima  are  shifted  in  altitude  and  smeared  out. 
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Figure  12.  Comparison  of  standard  deviations  (a)  with  the  sum  of  the  square  of  the  amplitudes  of  all 
spectral  bins  calculated  by  the  Fourier  transformation  (b).  Panel  (c)  gives  the  ratio  of  Panel  (a)  divided  by 
Panel  (b).  For  details,  see  text. 


Thus,  at  the  lowermost  heights  of  ^25  km  the  ratio  shown 
in  Figure  12b  approaches  zero.  One  should  keep  this 
weakness  of  the  FT  in  mind  when  comparing  these  results 
with  those  of  the  MEM/HA. 

[51]  To  give  a  realistic  impression  of  the  relation  of  GW 
variance  contained  in  the  several  wavelength  ranges,  visi¬ 
bility  corrections  as  deduced  in  section  2.4  are  applied.  The 
FT  wavelength  bins  are  corrected  with  scaling  factors 
according  to  their  center  wavelengths,  the  amplitudes 
obtained  from  MEM/HA  are  scaled  individually  as 
described  in  section  3.2.  In  contrast  to  the  FT  the  wave¬ 
length  ranges  of  zonal  means  calculated  from  MEM/HA 
results  can  be  chosen  arbitrarily.  In  Figure  13  we  compare 
the  results  of  the  FT  (uppermost  row)  to  those  obtained  by  a 
MEM/HA  with  interval  lengths  of  9  points  (13  km;  middle 
row)  and  5  points  (7.5  km;  lowermost  row).  The  left  column 
compares  the  FT  wavelength  bin  centered  at  20  km  (Panel 
(a))  with  MEM/FIA  results  in  the  range  from  15  km  to  30  km 
vertical  wavelength  (Panels  (d)  and  (g)).  The  middle  and 
right  columns  compare  the  sum  of  all  FT  bins  with  center 
wavelengths  shorter  than  10  km  (b)  and  6.7  km  (c), 
respectively,  to  MEM/HA  results  in  the  wavelength  ranges 
of  3.0-15.0  km  (e,  h)  and  of  3.0-10.0  km  (f,  i).  Good 
agreement  between  FT  and  MEM/HA  variances  is  found  for 
altitudes  above  30  km.  However,  the  deviations  below  30 
km  are  as  expected  because  the  analysis  interval  of  the 
MEM/HA  is  much  smaller  than  the  sliding  window  of  the 
FT  and  the  vertical  smearing  effects  are  much  less  of  a 
problem  using  MEM/HA.  Particularly  low  smearing  effects 
are  found  (cf  standard  deviations  in  Figure  12)  when  we 
reduce  the  window  length  of  the  HA  from  9  points  (13  km; 
middle  row)  to  5  points  (7.5  km;  lowermost  row).  The 
comparison  shows  that  a  7.5  km  interval  length  results  in 
reliable  zonal  mean  values  even  in  the  15-30  km  vertical 
wavelength  range. 

[52]  It  should  be  noted  that  both  methods  exhibit  a 
pronounced  tropical  maximum  and  subtropical  minima  in 
the  lower  stratosphere  only  for  short  vertical  wavelengths. 
This  is  an  important  issue  when  comparing  GW  climatol¬ 


ogies  based  on  different  satellite  instruments  [Preusse  et  al. , 
2000],  The  fact  that  we  observe  the  tropical  maximum  only 
at  short  vertical  wavelengths  may  be  related  to  the  higher 
occurrence  of  low  frequency  and  long  horizontal  wave¬ 
length  waves  in  the  equatorial  region  [Alexander  et  al, 
2002].  These  long  horizontal  wavelengths  in  the  tropics 
induce  errors  when  the  scaling  function  is  applied.  The 
deduction  of  the  scaling  function  is  based  on  the  assumption 
that  the  horizontal  wavelength  of  the  waves  is  of  the  order 
of  1000  km  or  shorter.  This  is  not  generally  true  in  the 
tropics.  At  low  vertical  wavelengths  the  reduced  sensitivity 
is  mainly  due  to  the  retrieval  process.  Thus,  it  does  not 
depend  on  the  horizontal  wavelength,  as  can  be  seen  from 
Figures  4  and  6.  At  long  vertical  wavelength  the  radiative 
transfer  causes  reduced  sensitivity  values  and  a  strong 
dependence  on  the  horizontal  wavelength  is  observed.  This 
leads  to  an  overcorrection  of  the  waves  with  both  long 
horizontal  and  long  vertical  wavelengths,  and  hence  to  an 
overestimation  of  the  tropical  wave  activity  shown  in  the 
left  column  of  Figure  13. 

[53]  From  the  results  presented  in  Figures  12  and  13  we 
conclude  that  the  MEM/HA  provides  reliable  information 
on  GW  temperature  variance  for  variable  choices  of  wave¬ 
length  intervals.  It  is  superior  to  the  FT  method  in  dealing 
with  nonstationary  wavelengths  and  amplitudes  in  profiles. 
We  also  find  that  a  9-point  sliding  window  is  sufficient  for 
analyzing  waves  with  wavelengths  longer  than  20  km.  The 
reliability  of  the  results  is  further  supported  by  the  excellent 
agreement  between  results  from  GPS/MET  [Tsuda  et  al., 
2000]  and  MEM/HA  amplitudes  inferred  from  CRISTA 
data  [Preusse  et  al.,  2000]. 

4.  Case  Study:  Stratospheric  GWs  Measured 
Above  South  America 

[54]  Previous  investigations  of  CRISTA- 1  temperature 
data  by  Preusse  et  al.  [1999]  and  Eckermann  and  Preusse 
[1999]  have  shown  that  on  day  310  (6  November  1994) 
very  large  amplitude  GWs  were  observed  over  the  southern 
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Figure  13.  Comparison  of  zonal  mean  spectral  power  in  CRISTA  residual  temperatures  using  FT  (upper 
row)  and  MEM/HA  on  a  9  point  (middle  row)  and  5  point  (lowermost  row)  sliding  window  for  different 
wavelength  ranges  (different  columns).  In  general,  the  agreement  is  good,  but  the  MEM/HA  provides 
better  altitude  resolution  and  thus  reproduces  the  tropical  maximum  at  more  realistical  altitudes  (compare 
also  Figure  12;  for  details,  see  text). 


tip  of  South  America.  Eckermann  and  Preusse  [1999] 
attributed  the  enhanced  perturbations  to  MWs  forced  by 
flow  over  the  Andes.  They  compared  a  time  series  of  four 
consecutive  days  of  CRISTA- 1  data  from  this  region  to 
hindcasts  of  the  Naval  Research  Laboratory  Mountain  Wave 
Forecast  Model  (NRL-MWFM).  It  has  been  shown  that  the 
model  was  able  to  reproduce  the  day-to-day  variations  in 
GW  activity.  However,  the  NRL-MWFM  results  were  con¬ 
centrated  over  the  mountain  ridge,  while  enhanced  fluctua¬ 
tions  in  the  CRISTA  data  extend  further  to  the  lee  side  and 
beyond  the  very  southern  tip  of  South  America  toward  the 


ocean.  In  addition  there  was  a  study  by  Tan  and  Eckermann 
[2000]  based  on  a  two-dimensional  nonlinear  high  resolu¬ 
tion  model.  However,  to  address  the  full  horizontal  structure 
of  the  GW  field  in  this  region  a  different  kind  of  simulation 
has  to  be  used.  We  therefore  simulate  waves  above  South 
America  on  this  day  by  employing  a  three-dimensional 
mesoscale  model. 

4.1.  Comparison  With  the  MM5  Mesoscale  Model 

[55]  The  model  we  use  is  the  fifth  generation  Pennsylva¬ 
nia  State  University  (PSU)/National  Center  for  Atmospheric 
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Figure  14.  GW  activity  over  South  America  on  6  November  1994  at  a  pressure  surface  of  30  hPa.  The 
vertical  wavelength  measured  by  CRISTA  and  the  viewing  direction  of  the  instrument  are  shown  in  Panel 
(a).  The  wavelength  information  is  significant  only  for  waves  with  amplitudes  ^1.5  K.  The  wave 
amplitudes  are  shown  in  Panel  (b)  together  with  the  phase  of  the  waves,  which  is  indicated  by  the 
direction  of  the  arrows.  An  upward-pointing  arrow  indicates  0°  and  increasing  phase  is  indicated  by 
counterclockwise  rotation.  Panel  (c)  shows  the  temperature  field  at  30  hPa  as  simulated  by  MM5  and 
Panel  (d)  shows  MM5  results  using  a  nested  grid  to  increase  the  horizontal  resolution.  Only  the  ran  with 
enhanced  resolution  exhibits  pronounced  wave  structures  at  the  red  dots,  which  give  the  location  of  the 
three  strong  amplitude  profiles  visible  in  Panel  (b).  For  details,  see  text. 


Research  (NCAR)  mesoscale  model  MM5  [Dudhia,  1993; 
Grell  et  al,  1994].  The  model  integrates  the  fully  compres¬ 
sible  set  of  fluid  equations  in  a  rotating  frame  of  reference. 
In  the  nonhydro  static  version  used  here  velocity,  pressure 
perturbation  and  temperature  are  the  prognostic  variables. 
The  model  domain  is  centered  around  65°W  and  45°S  with 
a  horizontal  domain  size  of  3000  km  x  3000  km.  A  local 
grid  refinement  (nested  domain  of  1210  km  x  1210  km) 
covers  the  mountain  ridge  on  the  west  side  of  South 


America.  The  horizontal  mesh  size  of  the  outer  domain  is 
Ax  =  30  km,  the  grid  length  of  the  nested  domain  is  10  km. 

[56]  The  vertical  resolution  z\z  is  ^260  m.  In  total,  99 
levels  are  used  up  to  the  model  top  at  10  hPa  (~28  km), 
where  a  radiative  boundary  condition  avoids  spurious 
downward  reflection  of  vertically  propagating  GWs.  Radi¬ 
ative  and  moist  processes  are  switched  off  as  the  prime 
concern  is  the  dynamics  of  MWs  at  upper  levels.  The  initial 
conditions  and  boundary  values  of  the  model  integration 
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were  prescribed  by  6-hourly  analyses  of  the  European 
Center  for  Medium-Range  Weather  Forecasts  (ECMWF) 
with  a  horizontal  resolution  of  1.125°  in  latitude  and 
longitude  and  15  levels  between  the  surface  and  10  hPa. 
This  model  setup  has  been  successfully  used  for  simulating 
stratospheric  waves  above  Scandinavia  and  southern  Green¬ 
land  [Dornbrack  et  al,  1998,  1999;  Leutbecher,  1998; 
Dornbrack  and  Leutbecher,  2001],  The  model  simulation 
starts  on  5  November  at  0000  GMT.  Data  shown  here  are 
simulated  for  6  November  0600  GMT,  which  is  closest  to 
the  CRlSTA-observations.  Thus,  the  model  had  already  run 
for  30  h,  when  the  data  shown  here  were  extracted. 

[57]  Figure  14  compares  the  CRISTA  data  with  the 
results  from  the  MM5  model  results  at  a  pressure  level  of 
30  hPa  (^24  km  altitude).  Figure  14a  shows  a  map  of 
MEM/HA  temperature  perturbations  from  CRISTA.  The 
color  code  gives  the  vertical  wavelength,  and  the  alignment 
of  the  lines  indicates  the  viewing  direction  of  the  telescope. 
Figure  14b  shows  the  corresponding  amplitudes  of  the 
waves  and  the  arrow  directions  indicate  their  phases.  The 
phases  are  determined  by  the  MEM/HA  from  the  vertical 
wave  structure  with  an  upward-pointing  coordinate  and  z  = 
0  at  the  plotting  altitude.  A  northward-pointing  arrow 
indicates  0°  and  increasing  phase  is  indicated  by  counter¬ 
clockwise  rotation  of  the  arrow.  (This  means:  90°  =  west¬ 
ward-pointing,  180°  =  southward-pointing,  270°  = 
eastward-po  inting) . 

[58]  If  the  CRISTA  sampling  is  sufficiently  dense  in  the 
horizontal  to  resolve  the  direction  of  the  phase  turn,  we  can 
infer  the  tilt  of  the  wave  fronts  in  the  vertical,  because  these 
phase  estimates  contain  vertical  information.  For  example,  a 
phase  increasing  from  west  to  east  (counterclockwise  rota¬ 
tion)  indicates  a  westward  tilt  of  the  phase  fronts.  The 
underlying  black  crosses  indicate  the  locations  of  the 
tangent  points.  Panels  (c)  and  (d)  give  results  from  hindcast 
simulations  using  the  MM5  model  which  will  be  discussed 
below. 

[59]  On  the  upwind  (westward)  side  of  the  Andes 
negligible  amplitudes  are  observed  in  Figure  14b.  The 
corresponding  wavelength  information  in  Figure  14a  is 
insignificant  for  profiles  with  amplitude  T  <  1  K.  Enhanced 
GW  amplitudes  are  observed  above  and  downwind  (east¬ 
ward)  of  the  landmass,  but  not  further  northward  than 
^40°  S.  This  can  be  explained  by  considering  the  disper¬ 
sion  relation.  For  a  two-dimensional  hydrostatic  MW  in  a 
constant  background  atmosphere  and  neglecting  the  Cori¬ 
olis  force  the  theoretical  vertical  wavelength  is  given  from 
(12)  by 


Figure  15.  Contours  of  the  upper  limit  of  the  vertical 
wavelength  of  MWs  (labels  in  kilometers)  on  6  November 
1994  at  30  hPa.  The  values  are  calculated  from  UKMO 
winds  and  CRISTA  background  temperatures  for  the  region 
from  30°S-60°S  and  30°W-90°W.  Overplotted  are  wind 
vectors  indicating  the  wind  direction  and  wind  speed. 
Highest  wind  speeds  shown  (largest  arrows)  are  60  m  s~'. 


vertical  wavelengths  too  short  to  be  detected  by  CRISTA 
(see  Figure  6). 

[60]  Enhanced  GW  activity  is  observed  in  the  lee  of  the 
Andes  and  above  the  Falkland  Islands.  However,  the 
vertical  wavelengths  observed  there  are  ^20  km  and  there¬ 
fore  a  factor  ^2  larger  than  (14)  predicts.  This  indicates  that 
the  source  of  the  observed  waves  is  probably  nonoro- 
graphic.  The  observation  of  nonorographic  waves  near  an 
isolated  island  on  the  southern  hemisphere  suggests  paral¬ 
lels  to  the  work  of  Guest  et  al.  [2000],  who  described  wave 
generation  by  a  frontal  system  at  the  upwind  side  of 
Macquarie  Island  (55°S). 

[61]  The  strongest  GW  amplitudes  are  observed  at  around 
43°S.  There  are  three  consecutive  profiles  with  high  ampli¬ 
tudes  and  alternating  phase.  These  three  profiles  have  been 
discussed  in  detail  by  Eckermann  and  Preusse  [1999].  They 
observed  coherent  wave  oscillations  at  altitudes  between 
15  km  and  30  km  and  inferred  that  these  waves  encountered 
a  critical  level  at  ^43  km  and  break  strongly  at  altitudes  just 
below.  The  vertical  wavelength  as  well  as  the  amplitude 
growth  with  altitude  provided  strong  evidence  that  these 
waves  are  MWs.  A  monochromatic  sinusoidal  wave  can  be 
expressed  as  [Preusse  et  al,  2001b] 

T'(x,z,  t)  =  7’(z)sm(G:  +  mz  —  tot  +  (j))  (15) 


(2Tt)^  Aid’ll^ 

^  N^/lP-  -  l/(4//2)  ~  A72 


(14) 


with  N  the  buoyancy  frequency  and  H  the  scale  height.  If 
the  angle  $  between  the  wind  direction  and  the  horizontal 
wave  vector  differs  from  180°,  i.e.,  if  the  wind  direction 
is  not  perpendicular  to  the  wave  fronts,  \z,2Z)  has  to  be 
scaled  by  a  factor  |cos  $|  given  the  approximation  on  the 
right  hand  side  of  (14).  Figure  15  plots  the  maximum 
vertical  MW  wavelength  calculated  from  UKMO  analysis 
winds  on  this  day  using  (14).  The  arrows  show  the  wind 
velocities,  which  are  directed  nearly  zonally.  Figure  15 
shows  that  northward  of  ^40°S  possible  MWs  have 


where  T'  is  the  residual  temperature,  t  denotes  the 
temperature  amplitude,  k  -  2tt/\j-  is  the  horizontal  wave 
number,  m  =  is  the  vertical  wave  number,  and  w  = 

2tt/t  is  the  ground-based  frequency.  For  an  observer  at  the 
ground  u)  =  0  in  the  case  of  a  MW.  A  phase  variation  in  the 
vertical  between  two  profiles  sampling  the  same  wave 
therefore  indicates  the  same  phase  variation  in  the 
horizontal.  This  is  also  discussed  by  Eckermann  and 
Preusse  [1999]  and  Preusse  et  al.  [2001b]  in  some  detail. 
Thus,  from  Figure  14b  one  can  conclude  that  the  horizontal 
wavelength  along  the  orbit  track  is  either  =  400  km,  in 
which  case  the  phase  changes  by  180°  between  the 
subsequent  profiles,  or  X^:  =  130  km,  in  which  case  the 
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phase  changes  by  540°  between  the  subsequent  profiles. 
The  latter  would  mean  that  CRISTA  undersamples  the  wave 
and  it  seems  that  the  question  of  whether  =  130  km  or 
400  km  cannot  be  answered  on  basis  of  these  CRISTA  data 
alone. 

[62]  However,  this  question  is  directly  related  to  the 
sensitivity  discussion  of  section  2.  At  this  location  the  main 
ridge  of  the  Andes  follows  the  coast  line  and  the  wind 
direction  is  nearly  perpendicular  to  the  ridge  at  all  altitudes 
{Tan  and  Eckermann,  2000].  We  therefore  can  assume  that 
the  wave  crests  are  aligned  with  the  coast  line  as  well.  This 
assumption  is  supported  by  the  measured  vertical  wave¬ 
length,  which  is  in  excellent  agreement  with  the  two-dimen¬ 
sional  dispersion  relation  (14)  [Eckermann  and  Preusse, 
1999]  and  therefore  indicates  a  two-dimensional  geometry. 
From  Figure  14a  we  see  that  the  three  profdes  with  the  high 
amplitudes  are  measured  at  an  angle  of  about  45°  between 
the  line  of  sight  and  the  coast  line  and  hence  the  wave 
crests.  Only  a  few  kilometers  apart  are  measurements  from  a 
second  orbit  where  the  line  of  sight  was  nearly  perpendic¬ 
ular  to  the  mountain  ridge.  From  Figure  14b  we  see  that  the 
amplitudes  measured  during  this  orbit  are  vanishing  (<1  K). 
The  two  sets  of  measurements  were  taken  at  0624  GMT  and 
0453  GMT  respectively.  Using  the  same  approximations  as 
for  (14)  the  vertical  group  velocity  of  a  GW  is  given  by 
[e.g.,  Andrews  et  al,  1987] 


and  is  ^  1  km  h“ '  for  a  wave  of 400  km  horizontal  and  6  km 
vertical  wavelength.  Thus,  we  do  not  expect  the  wave  field 
to  alter  sufficiently  during  this  short  time  difference  between 
the  two  measurements  to  explain  these  large  differences  in 
the  observed  amplitudes.  Thus  we  can  only  explain  the 
different  measured  amplitudes  if  the  horizontal  wavelength 
is  just  at  the  edge  of  the  sensitivity  distribution  described  in 
Figure  6.  This  is  the  case  if  we  assume  a  horizontal 
wavelength  of  130  km,  which  is  not  visible  to  CRISTA  if 
viewed  perpendicular  to  the  wave  crests,  but  becomes 
visible  when  the  wavelength  along  the  line  of  sight  is 
elongated  to  more  than  200  km  by  viewing  the  wave  crests 
at  a  significant  angle  a  (see  section  2.5  and  Figure  9).  This 
interpretation  is  further  supported  by  investigations  with  a 
two-dimensional  high  resolution  model  [Tan  and  Ecker¬ 
mann,  2000]. 

[63]  Figures  14c  and  14d  show  the  MM5  temperature  field 
in  the  outer  domain  at  30  hPa  at  0600  GMT  on  6  November 
1994  from  two  different  simulations:  Figure  14c  without. 
Figure  14d  with  a  local  grid  refinement.  The  three  red  points 
give  the  locations  of  the  profiles  investigated  by  Eckermann 
and  Preusse  [1999].  Only  the  higher  resolution  simulation 
exhibits  significant  stratospheric  GW  activity  above  these 
three  points.  Viewed  optimistically.  Profiles  1  and  2  are  each 
located  in  Figure  14d  at  the  eastern  edge  of  a  cold  trough, 
indicated  by  deep  blue,  thus  separated  by  about  one  wave¬ 
length.  Profiles  2  and  3  are  separated  by  a  third  cold  trough 
and  seem  to  be  separated  by  one  and  a  half  wavelengths.  For 
this  interpretation  one  has  to  assume  that  the  splitting  of  the 
large  blue  shaded  trough  spanning  from  70°W,  40°S  to 
60°W,  50°S  results  in  separate  troughs.  If  we  do  so,  the 
phase  fronts  are  aligned  parallel  to  the  coast  lines  and  waves 
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Figure  16.  Enlargement  of  the  region  in  the  vicinity  of 
Tierra  del  Fuego.  Amplitude  and  phase  infonnation  (cf 
Figure  14b)  are  overplotted  on  MM5  high-resolution  data 
(cf  Figure  14d).  To  guide  the  eye,  MM5  phase  fronts  are 
extended  to  the  south  (red  lines).  The  phases  of  the  CRISTA 
profiles  indicate  that  these  turn  southward  toward  a 
longitudinal  direction  on  the  south  coast  (red  dashed  lines). 


detectable  by  CRISTA  emerge.  A  wave  front  following  the 
main  direction  of  the  large-scale  trough  results  in  too  short  a 
vertical  wavelength  and  is  not  compatible  with  the  CRISTA 
observations.  Thus  the  MM5  model  indicates  a  short  hori¬ 
zontal  wavelength  of  wHO  km,  too.  However,  the  ampli¬ 
tudes  in  the  MM5  results  seem  to  be  strongest  in  the  vicinity 
of  the  main  mountain  ridge  (westward  of  70°W)  and 
decreasing  amplitudes  are  observed  at  the  east  coast  (profile 
3).  This  is  much  less  pronounced  in  two-dimensional  model 
simulations  [Tan  and  Eckermann,  2000]  and  not  found  in  the 
CRISTA  measurements  of  the  three  profiles.  Probably  the 
model  resolution  will  have  to  be  further  increased  for  this 
region.  This,  however,  requires  a  higher  resolution  topo¬ 
graphic  data  set  within  the  MM5  model  than  the  one  which 
was  available  to  the  authors. 

[64]  The  question  of  whether  the  wavelength  of  the  GW 
observed  near  42°S  is  more  likely  to  be  130  km  rather  than 
400  km  is  not  purely  academic.  As  has  been  pointed  out  by 
Eckermann  and  Preusse  [1999],  this  single  wave  event 
carries  a  large  part  of  the  total  momentum  flux  for  this 
latitude  and  a  130  km  wavelength  wave  would  lead  to  three 
times  larger  momentum  flux  deposition  (for  a  given  temper¬ 
ature  amplitude)  than  a  400  km  wave. 

[65]  Good  agreement  between  CRISTA  data  and  the 
MM5  results  is  found  at  the  southernmost  tip  of  South 
America.  The  CRISTA  data  in  Figure  14b  show  large 
amplitudes  in  the  vicinity  of  Tierra  del  Fuego,  reaching 
down  to  the  southernmost  boundary  of  measurements  at 
57°S.  An  enlarged  map  of  this  region  is  shown  in  Figure  16. 
The  CRISTA  profiles  are  overplotted  on  the  high-resolution 
MM5  temperature  field.  In  the  CRISTA  data  we  observe  a 
phase  front  along  the  west  coast  with  phases  around  0° 
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theoretical  [km] 

Figure  17.  Scatterplot  of  measured  versus  theoretical 
vertical  wavelength  in  the  area  42°S-50°S.  Symbols  show 
data  from  different  days,  from  day  309  (5  November)  to  day 
315  (11  November).  Almost  all  the  data  points  lie  within  the 
error  bars  (dashed  lines).  The  shading  indicates  amplitude- 
weighted  point  density.  For  details,  see  text. 


(three  northward-pointing  arrows),  a  second  one  200  km  to 
the  east  with  phases  around  180°  (three  southward-pointing 
arrows)  and  a  third  front  again  200  km  further  to  the  east  at 
about  0°.  (Note  that  at  55°S  two  profiles  in  the  first  and 
second  phase  front  are  overplotted  on  top  of  two  other 
profiles  with  nearly  the  same  orientation  and  a  similar 
amplitude,  respectively.)  The  profile  at  ^50°S,  73°W  is 
not  consistent  with  a  stationary  MW  as  the  long  vertical 
wavelength  >  20  km  indicates. 

[66]  A  very  similar  structure  is  visible  in  the  MM5  results. 
We  find  a  cold  trough  eastward  of  the  coast  followed  by 
warm  and  cold  phases  separated  by  ^200  km.  The  hori¬ 
zontal  wavelengths  retrieved  from  the  measurement  and  the 
model  are  in  good  agreement.  To  make  comparison  easier 
the  three  red  lines  show  the  extension  of  the  MM5  phase 
fronts  to  the  south. 

[67]  The  three  CRISTA  profiles  at  ^67°W  are  aligned 
nearly  meridionally  with  even  a  slight  tilt  to  the  north  east 
direction.  The  phase  fronts  produced  by  the  MM5,  however, 
are  aligned  approximately  parallel  to  the  coast  line.  Thus,  if 
the  southernmost  profile  is  located  just  between  trough  and 
crest  the  northern  profile  is  located  in  the  trough.  The 
phases  observed  by  CRISTA  show  this  behavior.  The 
phases  change  from  zero  (southernmost  profile)  to  a  neg¬ 
ative  value  (northern  profile,  slightly  eastward  direction). 
Thus,  the  comparison  of  the  CRISTA  and  the  MM5  data 
shows  that  the  model  agrees  with  the  data  in  the  horizontal 
wavelength  (^400  km)  and  in  the  direction  of  the  phase 
fronts.  The  amplitudes  measured  by  CRISTA  in  this  region 
are  about  4-5  K.  This  compares  very  well  with  peak-to- 
peak  variations  of  ^10  K  in  the  MM5  results.  The  only 
difference  is  that  the  CRISTA  measurements  indicate  wave 
trains  propagating  long  distances  downstream  (eastward) 
over  the  open  ocean  down  to  57°S,  which  is  less  pro¬ 


nounced  in  the  MM5  model  results.  It  further  appears  that  at 
the  south  coast  of  Tierra  del  Fuego  the  phase  fronts  of  the 
measured  waves  turn  toward  a  longitudinal  direction  (see 
the  red  dashed  line  in  Panel  6,  cf  the  phase  direction  of  the 
three  profiles  on  the  west  coast). 

[68]  From  Figure  14a  we  see  that  the  vertical  wavelength 
of  the  waves  at  55°S  is  about  10  km.  The  wind  speeds 
shown  in  Figure  15  are  directed  slightly  southeastward  and 
thus  blow  at  an  angle  $  «  40°  to  the  wave  fronts  at  55°S. 
We  therefore  expect  the  measured  wavelength  to  be  a  factor 
of  cos(40°)  «  0.77  shorter  than  the  wavelength  predicted  in 
Figure  15,  which  indeed  we  find. 

[69]  Eckermann  and  Preusse  [1999]  evaluated  the  verti¬ 
cal  wavelength  on  a  statistical  basis  for  the  region  north¬ 
ward  of  50°S.  We  reproduce  this  diagram  in  Figure  17.  The 
measured  vertical  wavelength  is  plotted  versus  the  predicted 
upper  limit  wavelength  using  (14).  Only  wavelength  esti¬ 
mates  for  oscillations  with  amplitudes  larger  than  2  K  are 
shown.  The  shading  gives  the  amplitude-weighted  point 
density.  The  uncertainty  of  the  retrieved  vertical  wavelength 
is  taken  from  the  scatter  in  Figure  5  (section  2)  and  is 
plotted  with  dashed  lines.  Only  a  very  few  points  lie  outside 
this  error  boundary.  Note  also  the  high  point  density  at 
^8  km  vertical  wavelength  along  the  1:1  line.  This  agree¬ 
ment  between  theoretical  and  observed  wavelengths  is 
evidence  that  the  waves  here  are  stationary  MWs  and 
furthermore  that  the  wave  fronts  here  are  aligned  approx¬ 
imately  perpendicular  to  the  wind  vector. 

[70]  Data  southward  of  50°S  are  plotted  in  Figure  18. 
Here,  as  discussed  above,  we  do  not  expect  the  values  to  be 
distributed  tightly  around  the  1 : 1  line  due  to  the  nonzero  $ 
for  waves  here.  We  therefore  draw  the  solid  line  \z,measured  = 
Q.n\,2D  onto  Figure  18.  As  expected  we  find  the  majority 
of  the  waves  now  centered  around  the  solid  line  and  inside 
the  respective  error  lines  (dashed),  indicative  of  $  ^  40°. 


Figure  18.  As  for  Figure  17,  but  for  50°S-57°S.  The  solid 
line  indicates  the  expectations  for  wave  fronts  at  an  angle 
<&  =  40°  to  the  mean  background  wind.  Most  points  cluster 
around  the  solid  line  inside  the  error  limits  (dashed  lines). 
The  dashed-dotted  line  indicates  \2d  plus  the  estimated 
precision  uncertainty  of  CRISTA.  For  details,  see  text. 
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However,  some  waves  may  stem  from  singular  mountains 
and  should  in  this  case  have  shorter  wavelengths  than  the 
dashed-dotted  line.  Nonetheless,  these  wavelength  estimates 
are  entirely  consistent  with  the  tilted  wave  fronts  simulated 
by  MM5  (Figures  14c,  14d,  and  16). 

4.2.  Interpretation  Employing  the  GROGRAT  Model 
4.2.1.  Questions 

[71]  The  close  comparisons  between  CRISTA  and  MM5 
data  in  Figures  14  and  16  show  conclusively  that  CRISTA 
accurately  measured  and  characterized  stratospheric  MWs 
over  the  Andes  on  6  November  1994.  To  gain  some 
physical  insights  into  the  MWs  in  both  the  CRISTA  data 
and  MM5  model  results,  we  conducted  detailed  ray-tracing 
simulations.  In  particular,  we  wanted  to  answer  the 
following  questions. 

[72]  (1)  Why  do  MWs  with  a  horizontal  wavelength 
^130  km  appear  to  dominate  at  40-45°S,  rather  than  the 
longer  ^400  km  wavelength  noted  further  to  the  south?  The 
Andean  ridge  remains  quite  wide  here  and  might  be 
expected  to  force  a  longer  wave  mode  throughout  [e.g., 
Tan  and  Eckermann,  2000].  Nonetheless,  it  does  not  arise 
here  in  either  the  data  [see  also  Eckermann  and  Preusse, 
1999]  or  MM5  results.  Exactly  why  is  unclear. 

[73]  (2)  Why  does  a  ^400  km  horizontal  wavelength 
wave  dominate  at  ^50-55°S,  and  why  does  it  have  a 
noticeable  northwest -southeast  alignment  of  its  phase  lines 
(Figures  14  and  16)?  Two  possible  explanations  come  to 
mind.  First,  the  topography  here  “hooks,”  from  north- 
south  alignment  to  a  more  northwest- southeast  aligned 
coastline  at  the  southernmost  tip  of  South  America  (e.g., 
Figure  14  of  this  paper  and  Plate  1  of  Tan  and  Eckermann, 
[2000]).  This  might  be  expected  to  force  a  plane  MW  with 
northwest- southeast  aligned  phase  fronts,  which  then  prop¬ 
agates  directly  into  the  stratosphere. 

[74]  In  addition,  however.  Figure  15  shows  that  there  is 
strong  latitudinal  shear  in  the  stratospheric  zonal  winds  {dU! 
dy)  over  this  region.  The  ray-tracing  equation  [Dunkerton, 
1984;  Marks  and  Eckermann,  1995] 


reveals  that,  for  an  initially  westward  wave  vector  {k<0,l  = 
0),  the  strong  latitudinal  wind  shear  {dU/dy  <  0)  at  50-55°S 
in  Figure  15  can  lead  to  the  development  of  a  southward 
wave  number  component  (/  <  0)  through  lateral  refraction. 
This  yields  a  southwestward  wave  number  vector  {k<0, 1  < 
0)  and  thus  northwest- southeast  GW  phase  lines  orthogo¬ 
nal  to  this  wave  vector,  as  observed  in  Figures  14  and  16. 
Ray-tracing  simulations  allow  us  to  test  both  mechanisms  to 
some  extent. 

[75]  (3)  Why  do  the  longer  MWs  at  ^45-55°S  appear  to 
propagate  significant  distances  downstream  of  the  Andes? 
One  possible  reason  is  that  these  longer  wavelength  waves 
(small  k)  are  of  inertia  gravity  waveform,  such  that  \kU/f\  ^ 
0(1).  Such  low  frequency  MWs  can  leak  downstream  in 
two  and  three  dimensions  due  to  a  mismatch  between  their 
horizontal  phase  and  group  velocities  induced  by  Coriolis 
effects  [e.g.,  Dudhia,  1993].  In  three  dimensions,  however, 
even  hydrostatic  stationary  MWs  aligned  at  an  angle  to  the 
mean  flow  can  propagate  downstream  [e.g..  Smith,  1980]. 


Transience,  nonhydrostatic  effects  or  nonlinearity  can  also 
yield  downstream  penetration  of  wave  energy.  Ray-tracing 
can  be  used  to  test  the  significance  of  inertial  modifications 
and  three-dimensionality  to  downstream  penetration  of 
wave  energy  here. 

4.2.2.  Experiments 

[76]  Ray-tracing  was  performed  with  the  Gravity  Wave 
Regional  or  Global  Ray  Tracer  (GROGRAT)  Model  [Marks 
and  Eckermann,  1995;  Eckermann  and  Marks,  1997].  The 
GROGRAT  simulations  reported  here  included  refraction 
due  to  atmospheric  gradients  in  all  three  spatial  directions, 
but  did  not  include  time  dependence  of  the  background 
atmosphere.  Wave  amplitudes  were  solved  using  constancy 
of  vertical  wave  action  flux,  and  dissipation  due  to  infrared 
damping,  turbulent  damping  and  saturation  were  all 
included:  see  the  work  of  Marks  and  Eckermann  [1995] 
for  more  details.  The  ray-tracing  was  performed  through 
background  atmospheres  provided  by  analyses  from  NASA’s 
Data  Assimilation  Office  (DAO)  [Coy  and  Swinbank, 
1997]  and  the  National  Centers  for  Environmental  Pre¬ 
diction  (NCEP)  reanalyses  [Kalnay  et  al,  1996;  Kistler 
et  al.,  2001].  NCEP  and  DAO  ray -tracing  results  were 
broadly  similar:  we  present  the  DAO-based  results  here 
due  to  their  more  complete  stratosphere  and  some  questions 
about  the  reliability  of  NCEP  stratospheric  data  over 
the  Andes  [Trenberth  and  Stepaniak,  2002].  Results  from 
4  DAO  ray-tracing  experiments  for  6  November  1994  at 
0000  GMT  are  shown  in  Figure  19.  These  experiments  were 
devised  to  cast  light  on  the  questions  posed  above. 

[77]  In  all  4  cases,  stationary  MWs  were  launched  from 
locations  shown  with  green  diamonds  in  Figure  19  at  an 
altitude  of  2  km  (i.e.,  over  and  just  above  the  Andes)  with 
a  starting  horizontal  velocity  amplitude  of  1  m  s“^  Note 
that  no  attempt  was  made  here  to  assign  different  ampli¬ 
tudes  based  on  the  underlying  topography,  in  contrast  to 
MM5  results  and  the  MW  model  simulations  of  Ecker¬ 
mann  and  Preusse  [1999].  Waves  were  traced  to  an 
altitude  of  24  km,  which  is  ^30  hPa.  Final  amplitudes 
at  24  km  in  Figure  19  are  given  in  terms  of  the  nonnalized 
temperature  perturbation,  for  easier  comparisons  with 
Figures  14  and  16.  We  traced  two  types  of  waves,  one 
with  an  initial  horizontal  wavelength  of  130  km,  the  other 
with  a  400  km  wavelength,  results  for  which  are  shown  in 
the  left  and  right  columns  of  Figure  19,  respectively.  In 
both  cases  we  also  launched  the  waves  with  two  different 
azimuth  groupings.  In  the  first  case,  wave  vectors  were 
initially  aligned  westward  ((j))„„Y  =  0°.  These  results  are 
shown  in  the  top  row  of  Figure  19,  and  allow  us  to  assess 
the  role  of  latitudinal  refraction  in  setting  observed  phase- 
tilts  in  the  wave  patterns.  In  the  second  case,  southwest- 
ward  wave  vectors  were  imposed  on  the  waves  at  their 
launch  heights  (((j))„„Y  =  20°  and  40°),  to  mimic  MWs 
launched  by  two-dimensional  topography  with  a  similar 
tilted  alignment.  Results  for  these  waves  are  given  in  the 
bottom  row  of  Figure  19.  In  all  cases,  waves  with  vertical 
wavelengths  <5  km  at  the  final  height  of  24  km  were 
removed,  to  account  for  the  insensitivity  of  CRISTA  to 
these  small  vertical  scales. 

4.2.3.  Results 

[78]  The  ray-tracing  results  in  Figure  19  show  several 
features  that  resemble  both  the  CRISTA  data  and  MM5 
model  results  from  Figures  14  and  16. 
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(a)(  Uni,  =1 33km, (  (l)L,  =0  ° 


(c)(  Uini.  =400km,(  (I)),,  =0  ° 


Figure  19.  GROGRAT  simulations  of  stationary  MW  propagation  from  2  to  24  km  using  DAO 
“MAESAE”  assimilations  for  6  November  1994  at  0000  GMT.  See  text  for  further  details  on  these 
simulations.  Green  diamonds  show  starting  locations  at  2  km.  Thick  black  curves  show  the  group 
propagation  paths  of  those  waves  that  propagate  to  24  km  and  have  vertical  wavelengths  >5  km.  Vectors 
show  wave  number  direction  K/,  —  {k,  1)  at  24  km,  with  length  equal  to  total  normalized  peak 
temperature  amplitude  (see  also  color  code  and  color  bar).  Diameters  of  fdled  circles  are  proportional  to 
the  final  horizontal  wavelength  at  24  km.  Results  are  shown  for  MWs  of  (a)  initial  horizontal 
wavelength  =133  km  and  initial  zonal  (westward)  wave  number  direction  4) =  0°,  (b)  {\h)imt  = 

133  km  and  =  20°  and  40°  (southwestward),  (c)  (\h)i,nt  =  400  km  and  cj),-,,,-,  =  0°,  and  (d)  (\h)imt  ^  400 
km  and  =  20°  and  40°. 


[79]  First,  Figure  19a  shows  that  a  130  km  wavelength 
MW  can  reach  the  stratosphere  at  40-43°S,  whereas  400  km 
MWs  initially  aligned  either  westward  (Figure  19c)  or 
southwestward  (Figure  19d)  do  not  reach  the  stratosphere 
at  ~40-43°S.  This  occurs  in  the  GROGRAT  simulations 
because  low-level  winds  were  sufficiently  weak  that  wave 
frequencies  near  the  ground  were  less  than  the  local  inertial 
frequency  |/|:  i.e.,  \kU\  <  \f\  for  the  400  km  waves.  The 


130  km  waves,  however,  could  propagate,  since  their  larger 
k  values  satisfied  \kU\  >  \f\.  This  explains  why  a  400  km 
wave  did  not  emerge  in  the  stratosphere  at  these  latitudes, 
despite  the  presence  of  very  wide  mountain  ridges  at  ^40- 
50°S  that  could  plausibly  force  such  a  wave. 

[so]  Figure  19a  also  shows  that  ray  paths  of  the  130  km 
westward-aligned  MWs  barely  deviate  from  their  initial 
positions:  i.e.,  they  propagate  almost  purely  vertically.  This 
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indicates  that  the  waves  are  hydrostatic  and  not  affected  by 
either  rotational  effects  or  three-dimensionality  to  any  great 
degree  (note  that  imposition  of  three-dimensionality  via 
wave  number  tilts  in  Figure  19b  did  yield  some  downstream 
penetration,  as  expected  [e.g.,  Smith,  1980]).  Conversely,  the 
400  km  westward-aligned  wave  groups  in  Figure  19c  did 
propagate  significant  distances  to  the  south  and  east.  This  is 
due  to  two  effects.  The  first  is  inertial  modification  due  to 
intrinsic  frequencies  that  are  not  large  compared  to  |/|. 
Second,  a  southward  drift  occurs  due  to  development  of 
meridional  wave  number  components  /  due  to  lateral  refrac¬ 
tion  via  (17),  as  can  be  seen  from  the  tilted  final  wave 
number  vector  directions  plotted  in  Figure  19c. 

[81]  While  the  tilts  acquired  due  to  refraction  in 
Figure  19c  are  significant,  they  are  not  as  large  as  those 
seen  at  50-55°S  in  the  MM5  results  of  Figures  14  and  16. 
Thus,  some  initial  tilting  due  to  curving  topography  here 
also  appears  to  be  needed.  Figure  19d  shows  results  for  the 
400  km  waves  with  initial  tilts  of  20°  and  40°  imposed.  The 
results  here  most  closely  resemble  the  MM5  results  and 
CRISTA  data  at  50-55°S  in  Figures  14  and  16,  in  that  this 
long-wavelength  wave  energy  is  confined  mostly  to  lati¬ 
tudes  ^50-55°S,  has  significant  meridional  wave  numbers 
(and  hence  phase  tilts),  and  exhibits  greater  penetration 
downstream  of  the  Andes.  Interestingly,  it  produces  a  region 
of  very  significant  downstream  penetration  at  ^52-55°S 
that  resembles  that  seen  in  the  CRISTA  and  MM5  data 
(Figures  14  and  16).  Furthermore  the  waves  at  the  south 
coast  of  Tierra  del  Fuego  in  Figure  19d  have  a  much  smaller 
meridional  tilt  than  the  waves  slightly  northward.  Tierra  del 
Fuego  is  already  inside  the  polar  vortex  and  so  the  latitu¬ 
dinal  gradient  of  the  zonal  wind  is  much  smaller  or  even 
reversed  (cf.  Figure  15).  This  meridional  alignment  is  in 
agreement  with  the  CRISTA  results,  where  we  found  some 
indication  for  a  turn  to  meridional  alignment  (Figure  16,  red 
dashed  lines). 

[82]  In  suimnary,  then,  these  GROGRAT  simulations 
reveal  that  shorter  wavelength  waves  at  40-45°S  occur 
because  400  km  waves  cannot  propagate  from  the  Andes 
at  these  latitudes  due  to  weak  low-level  winds,  whereas 
shorter  130  km  waves  can  propagate  into  the  stratosphere. 
Further  south,  where  winds  are  stronger  throughout,  the 
400  km  wave  can  propagate  from  topography  and  emerge 
strongly  in  the  stratosphere,  as  observed.  Its  tilted  phase 
alignment  and  strong  downstream  propagation  appears  to 
result  mainly  from  the  curved  topography  in  this  south¬ 
ernmost  region  of  South  America,  but  it  is  also  accentu¬ 
ated  by  lateral  refraction  effects  via  (17)  and  inertial 
modifications. 

5.  Summary  and  Conclusions 

[83]  In  the  recent  past  a  large  number  of  GW  observations 
from  satellite  instruments  have  been  reported  and  upcoming 
missions  will  also  search  for  GW  signatures.  However,  so 
far  the  sensitivity  to  GWs  has  only  been  estimated  in  detail 
for  instruments  measuring  saturated  radiances  [Wu  and 
Waters,  1997;  McLandress  et  al.,  2000].  In  this  paper  a 
detailed  estimation  of  the  sensitivity  to  GWs  of  emission 
limb-sounding  under  optically  thin  conditions  is  presented. 
We  have  shown  that  one  needs  to  consider  the  radiative 
transfer  throughout  a  GW  field  in  two  dimensions.  An 


analytical  approximation  derived  in  section  2.2  enabled  us 
to  separately  discuss  the  limitations  due  to  the  radiative 
transfer,  which  are  universal  to  all  instruments  using  the 
same  technique,  and  the  limitations  due  to  the  sampling  and 
retrieval,  which  are  characteristic  of  the  CRISTA  instru¬ 
ment.  In  addition,  numerical  simulations  based  on  a  fast 
radiative  transfer  scheme  \Gordley  et  al.,  1994;  Marshall  et 
al,  1994]  allow  us  to  take  into  account  variations  of  the 
background  atmosphere  as  well  as  wave  amplitude  changes 
with  altitude. 

[84]  Due  to  the  radiative  transfer,  only  waves  with 
horizontal  wavelengths  longer  than  ^200  km  along  the 
limb  ray  can  be  detected.  An  unexpected  decrease  of 
sensitivity  to  short  horizontal  and  long  vertical  wavelengths 
is  found,  which  is  very  important  for  the  intercomparison  of 
recent  GW  climatologies  obtained  from  different  satellite 
instruments  {Alexander,  1998;  Preusse  et  al.,  2000].  The 
vertical  sampling  and  the  retrieval  suppresses  waves  with 
vertical  wavelengths  shorter  than  about  5  km. 

[ss]  The  sensitivity  varies  only  moderately  with  latitude 
and  altitude.  Thus,  geographical  distributions  can  be 
retrieved  accurately  from  temperature  data.  We  presented 
results  of  CRISTA’s  sensitivity  to  GWs  using  two  numer¬ 
ical  simulations:  one  with  and  one  without  amplitude 
growth  dependent  on  altitude.  Both  runs  agree  well  with 
each  other  and  with  the  results  from  a  simplified  analytical 
approach.  This  indicates  the  basic  generality  of  the  obtained 
sensitivity  estimates,  which  allows  us  to  apply  these  results 
to  other  infrared  limb  scanning  instruments.  Two  examples 
are  LIMS  {Fetzer,  1990;  Fetzer  and  Gille,  1994,  1996], 
which  has  been  used  for  previous  studies,  and,  in  future, 
HIRDLS,  which  will  provide  a  long  time  series  of  high 
spatial  resolution  data  and  which  therefore  ought  to  be  very 
well  suited  for  GW  studies  based  on  our  results  here.  For 
these  instruments  section  2  of  this  paper  presents  the 
sensitivity  estimations  imperative  for  any  satellite  GW 
analyses. 

[86]  In  addition,  for  horizontal  wavelengths  below  ^  1 000 
km,  the  sensitivity  does  not  depend  to  any  significant  extent 
on  the  degree  of  amplitude  growth  with  altitude.  For  these 
shorter  horizontal  wavelengths  (^300-800  km)  the  sensi¬ 
tivity  depends  almost  entirely  on  the  vertical  wavelength 
(Figure  6).  This  universal  behavior  of  the  sensitivity  allows 
us  to  deduce  a  one-dimensional  correction  factor,  which 
compensates  for  the  sensitivity  changes  as  a  function  of 
vertical  wavelength  [Preusse  et  al,  2000].  However,  the 
detection  limits  of  ^200  km  horizontal  wavelength  along 
the  LOS  and  ^5  km  vertical  wavelength  still  remain  and 
cannot,  of  course,  be  compensated  for. 

[87]  GW  signals  were  extracted  from  CRISTA  temper¬ 
atures  with  the  help  of  a  Kalman  filter  and  were  analyzed 
using  a  combination  of  MEM  and  HA.  This  MEM/HA 
method  compares  well  with  the  results  obtained  from  Four¬ 
ier  Transforms,  but  provides  a  better  altitude  resolution. 
Both  methods  show  that  a  low  altitude  tropical  maximum  in 
GW  variances  is  observed  in  short  vertical-wavelength 
GWs  only.  As  has  been  shown  recently  by  Preusse  et  al. 
[2000],  this  finding  explains  differences  in  the  patterns 
observed  by  different  satellite  instruments.  In  particular 
the  CRISTA  data  for  wavelengths  shorter  than  <  10  km 
agree  very  well  with  results  from  GPS/MET  [Tsuda  et  al, 
2000].  The  fact  that  we  observe  the  tropical  maximum  only 
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at  short  vertical  wavelengths  is  probably  due  to  a  relatively 
higher  occurrence  of  low  frequency  and  long  horizontal 
wavelength  waves  in  the  equatorial  region  [Alexander  et  al. , 
2002]. 

[88]  In  addition,  the  MEM/HA  analysis  provides  the 
vertical  wavelength  of  monochromatic  wave  trains  as  a 
function  of  altitude.  The  precision  of  these  retrieved  verti¬ 
cal-wavelength  determinations  is  estimated  to  be  approx¬ 
imately  25%  (see  Figure  5).  This  allows  us  to  compare  the 
wavelength  of  measured  waves  with  theoretical  predictions. 

[89]  The  high  spatial  resolution  of  the  CRISTA  measure¬ 
ments  is  utilized  in  a  case  study  of  MWs  observed  in  the 
stratosphere  over  the  southern  tip  of  South  America.  At 
about  43°S  three  profiles  with  very  high  amplitudes  were 
observed  on  6  November  1994.  With  the  help  of  the 
sensitivity  arguments  discussed  earlier  we  find  evidence 
that  the  observed  signatures  stem  from  a  MW  with  a 
horizontal  wavelength  of  ^130  km.  This  finding  is  sup¬ 
ported  by  numerical  modeling  with  a  two-dimensional  high- 
resolution  mesoscale  model  [Tan  and  Eckermann,  2000] 
and,  in  addition,  by  simulations  employing  the  GROGRAT 
ray-tracing  model  [Marks  and  Eckermann,  1995;  Ecker¬ 
mann  and  Marks,  1997].  The  observed  vertical  wavelengths 
are  for  all  observations  northward  of  50°N  in  good  agree¬ 
ment  with  the  theoretical  expectations  for  a  two-dimen¬ 
sional  geometry.  The  results  emphasize  that  a  detailed 
knowledge  of  the  radiative  transfer  is  necessary  to  fully 
understand  the  satellite  observations. 

[90]  At  the  southernmost  tip  of  South  America  we  find  a 
rich  field  of  waves  of  ^400  km  horizontal  wavelength. 
Here,  the  wave  crests  and  troughs  are  no  longer  oriented 
orthogonal  to  the  wind  field.  Thus,  the  vertical  wavelength 
is  shorter  than  predicted  in  a  two-dimensional  geometry  (cf. 
Figure  18).  The  CRISTA  observations  agree  with  the  MM5 
results  in  horizontal  wavelength  (^400  km),  amplitude 
(^5  K)  and  the  alignment  of  the  wave  fronts.  Ray-tracing 
calculations  with  the  GROGRAT  model  indicate  that  the 
orientation  of  the  wave  fronts  at  an  angle  of  ^40°  with 
respect  to  the  background  wind  is  predominantly  due  to  the 
hooklike  topography  at  the  southernmost  tip  of  South 
America,  although  latitudinal  shear  in  zonal  wind  is  also 
an  important  effect. 

[91  ]  The  observed  waves  extend  further  to  the  south  than 
in  the  MM5  simulations.  Still  the  agreement  between  the 
MM5  results  and  the  measurements  is  excellent  and  leaves 
little  doubt  that  CRISTA  is  accurately  measuring  and 
characterizing  mesoscale  MWs  in  this  region  of  the  strato¬ 
sphere.  It  should  be  noted  that  CRISTA  passes  this  region 
on  three  subsequent  orbits  and  so  the  wave  was  viewed  at 
different  viewing  angles.  The  results  from  these  different 
overpasses  agree  well  with  each  other,  showing  that  a 
horizontal  wavelength  of  \z  =  400  km  is  well  above  the 
sensitivity  limit  as  was  inferred  in  our  sensitivity  study  in 
section  2. 

[92]  The  horizontal  sampling  on  6  November  was  200 
km  along  track  and  600  km  across  track  for  a  single  pass  of 
the  satellite  [Riese  et  al.,  1999].  In  the  geographic  region  of 
the  case  study  the  horizontal  resolution  was  further 
increased,  because  this  region  was  located  at  the  southern 
turning  point  of  the  orbit  track.  Therefore,  the  same  region 
was  sampled  during  subsequent  orbits.  The  case  study 
demonstrates  that  this  very  high  horizontal  resolution  of 


the  measurements  is  essential  for  revealing  the  spatial 
structure  of  the  waves.  In  particular,  the  information  about 
the  three-dimensional  spatial  structure  provides  valuable 
test  cases  for  numerical  models.  The  comparison  of  the 
CRISTA  results  with  the  MM5  simulations  demonstrates 
that  the  model  captures  the  salient  features  of  the  GW  field 
when  run  at  a  high  resolution.  However,  the  case  study  also 
provides  evidence  that  a  200  km  sampling  along  track  is 
still  too  coarse  to  fully  resolve  the  horizontal  wave  struc¬ 
tures  detected  by  the  triple  telescope  CRISTA  instrument. 
The  detailed  investigations  of  small-scale  structures  by 
CRISTA  in  the  case  study  was  only  possible  due  to 
sampling  the  south  tip  of  South  America  in  several  sub¬ 
sequent  overpasses  and  is  therefore  a  result  of  using  three 
telescopes  simultaneously.  Thus  a  case  study  like  the  one 
presented  here  is  at  its  present  stage  possible  only  with  the 
triple-telescope  CRISTA  instrument  and  will  be  restricted  in 
the  near  future  to  those  instruments  providing  data  of  very 
high  spatial  resolution,  i.e.,  the  HIRDLS  instrument. 
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